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?t : INTRODUCTION 
> . 

^ , Post-diction has become an art form at AGS and SPS energies. The timing of QM99, 
^ [ just as the RHIG facility comes on-line after nearly two decades of planning and con- 

■ " " ' struction, provides a unique opportunity to break out of that mold. The one day Quark 
Matter 1999 session, organized by M. Gyulassy, had the goal of updating and document- 
ing the Pi?i?-dictions for the upcoming RHIG experiments. Of course, sticking one's neck 
out just as the experimental "guillotine" starts its descent is risky. The brave partici- 
pants in this session presented a wide range of predictions summarized in this collection 
of contributions. 

The contributions are organized into five section: (1) Global and Inclusive Observables 
that probe entropy production and initial conditions in AA, (2) Hadron Flavor Observ- 
ables that probe hadro-chemical equilibration and baryon number transport, (3) HBT 
and Gollective Flow Observables that probe the space-time volume of the freeze-out sur- 
face of the reaction as well as the equation of state of ultra-dense matter, (4) Jets and 
Penetrating Rare Probes which are sensitive to the highest density phase of the reaction, 
and (5) Exotic Possibilities that illustrate some of the more speculative novel phenomena 
that may be explored at RHIG. 



1. GLOBAL AND INCLUSIVE OBSERVABLES 
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1.1. K.J. Eskola: Lower limits for (Et) and {Net) from pQCD &z hydrodynamics 
at the central rapidity unit in central Au-Au collisions at RHIC [|iy2| 

In heavy ion collisions at very high cms-energies the initial, very early particle and 
transverse energy production at central rapidities is expected to be dominated by multi- 
ple production of minijets, i.e. gluons, quarks and antiquarks with transverse momenta 
Pt ~ 1...2GeV ^ Aqcd 0- Assuming independent multiple semi-hard parton-parton 
collisions, the average transverse energy carried by the minijets produced with > Po 
the central rapidity window Ay = 1 in a central (b = 0) AA-coUision can be computed 
in the lowest order (LO) perturbative QCD (pQCD) as [H 

^t^qcd(v^,Pt > Po, Ay, 0) = Taa{0) f dpTdyidy2- — Pt- (1) 

JpoAy dpTdyidy2 

The differential cross section above is that of each binary (LO) parton collision, 

, J = ^ ^lfi/A{Xl,Q)x2fj/A{x2,Q)^ 

dp^dyidy2 dt 

1,1,9 

where the rapidities of the outgoing partons are yi and 1/2 • The nuclear collision geometry 
is accounted for by the standard nuclear overlap function TaaW ~ A'^/7iR\. In order 
to have at least one of the two outgoing partons in each minijet collision in the rapidity 
window Ay, appropriate kinematical cuts have to be made, see for details. We use 
the GRV94-LO parton distributions 0, and, to consistently include both the x- and the 
Q-dependence of the nuclear effects in the parton distributions ( fi/A{x, Q) 7^ fi/p{x, Q) ), 
we use the recent EKS98-parametrization We emphasize that the results presented 
here are obtained by using merely the lowest order pQCD but in order to approximate 
the effects of the next-to-leading-order terms in the minijet cross sections, we include 
i^'-factors 2.0 (1.5) for RHIC (LHC), together with the scale choice Q = px- 

From saturation arguments |^ , and from requiring agreement with inclusive pion spec- 
tra at central rapidities in pp collisions at a/s = 200 GeV, we expect po ~ 1...2 GeV for 
Au-Au collisions at RHIC. The few-GeV minijets are produced within a short proper time 
Tj ~ l/po ~ 0.2. ..0.1 fm/c, so they serve as initial conditions for further evolution of the 
early QGP. In addition to the pQCD component, at RHIC we also expect a non-negligible 
non-perturbative (soft) component in the initial transverse energy production. For the 
simple estimates here, we take the soft component directly from the measured average 
at 7] ~ at central Pb-Pb collisions at the SPS |^. The initial transverse energy at 
T = 1/pq thus becomes E\~ = -Et,pQcd + E^^"^^, as shown by the dashed lines in Fig. la. 
Dividing by the initial volume Vi = irR^TiAy, we get a Bj or ken-estimate of the initial 
energy density: = E^/Vi. In a fully thermalized, 1+1 dimensional boost-invariant 
hydrodynamic system, the pdV work causes the energy density to decrease as e ~ r~^/^ 
0] and, especially, the transverse energy to decrease as E^ ~ r~^/^. We do not attempt 
to follow the system through a phase transition here but simply decouple the system at 
6/ = 1.5 GeVfm"'^. The resulting E^'^'^^ represents a lower limit of {E^) (mod the decrease 
in the mixed/hadronic phase), and is plotted by the solid curves in Fig. la. It should also 
be noted that when transverse expansion is included, the loss of i?T is less than in the 
1-1-1 dimensional case considered above. 
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Figure 1. Final {Et) and {Nch) in the central rapidity unit. 



To get a corresponding estimate of the final state charged pion multiplicity (Fig. lb), we 
convert the initial energy density into a temperature Tj, from which the initial rapidity 
density of entropy Si can be computed. For simplicity (in spite of gluon dominance), let 
us assume full thermalization of gluons and quarks here. In an isentropic boost-invariant 
1+1 dimensional fiow, the rapidity density of entropy is conserved, so Si = Sf ^ 4|Ai'ch 
0, where N^h is the final state charged pion multiplicity in the central rapidity unit. In 
other words, we obtain a lower limit of A^^ch from entropy of the initial state. 

To show the dependence on the transverse momentum cut-off, we plot the lower limits 
for (-Et) and N^h as functions of po (the solid curves). The upper limit of E^^^^ is the 
initial E^, so po can be fixed from some other arguments, the measured (E^) serves in 
principle also as a measure of thermalization. Our favorite estimates for the lower limits 
can be read off from the figures at pq ~ 1.5 for RHIC and at po ~ 2 GeV for the LHC. 
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1.2. Yu.V. Kovchegov: Nuclei and Classical Yang Mills 

This contribution is based on tlie paper in collaboration with A. H. Mueller. The 
transparencies can be viewed at the web address given in Q]. 

First let us consider scattering of a gauge invariant current j = ~\F^j^Fl^i^ on a nucleus. 
The single gluon inclusive production amplitude for this process viewed in covariant gauge 
d ■ A = is depicted in Fig. ^ The current (dashed line) interacts with a nucleon in the 
nucleus producing a gluon, which then rescatters on the other nucleons in the nucleus. 
The interactions of the propagating gluon with nucleons are pictured in the quasi-classical 
approximation That limits the interactions to one (inelastic) or two (elastic) gluon 

exchanges (see Fig. ^ in the eikonal approximation. Interactions involving more gluons 
are suppressed by the powers of a^. We consider the case of not very high center of mass 
energy and not very large transverse momentum of the produced particles, so that there 
is no large logarithms to enhance those extra powers of as- Our expansion is equivalent 
to resummation in the powers of the effective parameter a'^A^^'^ HJl]; with A the atomic 
number of the nucleus. 

A calculation outlined in [|l| yields for the number of gluons produced at the transverse 
position x±: 



where we assumed that the nucleus is a sphere of radius R and has a constant density p. 
The gluon distribution of each nucleon is given by xG{x, Q^) and is taken at the two-gluon 
level [jl| . We note that in covariant gauge everything is given by the final state interactions, 
which result in a simple Glauber expression of Eq. (|^). If one views the process in the 
light cone gauge A^ = 0, where the nucleus is moving in the "plus" direction, everything 
can be pictured in terms of initial state interactions. The nuclear wave function, given by 
non-Abelian Weizsacker- Williams field , provides us with a gluon which already has all 
the information about the multiple rescatterings that broaden its transverse momentum. 
This gluon simply interacts with the current producing a final state on-shell gluon [Q]. 
The resulting production cross section is, of course, the same as given by Eq. (|^), since 
it is a gauge invariant object. 

Similar effects happen in the proton-nucleus (pA) collisions. We are just going to state 
the answer for the total inclusive gluon production cross section in the quasi-classical 



exp 



m - 1 



pxG{x, 



(3) 




Figure 2. Gluon production in the current-nucleus interaction as envisioned in the text. 
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approximation 
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Here (l^) is the typical transverse momentum squared of the produced gluon For- 
mula (^) has a very simple physical interpretation: the first term on the right hand side 
corresponds to the gluon emission off the proton after the proton passes through the 
nucleus and the second term provides us with gluon's transverse momentum broadening 
corresponding to the case when the gluon is already present in the proton's wave func- 
tion before the collision and multiply rescatters on the nucleons in the nucleus during the 
collision. 

Eq. (^) is written for the case of <^ {^D- The general formula, which is given 
in |T| is plotted in Fig. ^ for the case of proton-gold collision, where we assumed a 
simple superposition form of the gluon distribution function of a nucleon xGp^n ~ SxGq = 
ln((5^/A^). In Fig. |^ we show the logarithm of the gluon production cross section 
as a function of the transverse momentum of the gluon /. We put (Z^) = 2 GeV^ for 
RHIC is, in general, energy dependent). At this lowest order the rest of Eq. (H) is 
energy (and rapidity) independent. 

It is a very interesting and challenging problem to calculate the gluon production cross 
section for the case of nucleus-nucleus collision in the quasi-classical approximation, i.e., 
resumming all powers of a'^A^^^ (strong field limit). Some attempts in that direction have 
been made already [Q. However they correspond to the weak field limit IHJ^. The full 
nucleus-nucleus problem is still to be solved. 
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Figure 3. Differential gluon production cross section as a function of gluon's transverse 
momentum for p+Au collision. 
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1.3. X.N. Wang: HIJING at RHIC 

In this section we will present estimates of hadron spectra in central Au+Au collisions at 
RHIC energy from perturbative QCD parton models (for high px spectra) and HIJING 
model (for dNc\i/dy). Since there are many nuclear effects (including the formation of 
QGP) which we still don't quantitatively understand, the extrapolation from pp data 
at the RHIC energy and p + A, A + A data at SPS energy is inevitably rigged with 
uncertainties. We will try to estimate these uncertainties by considering different scenarios 
of various nuclear effects, like nuclear parton shadowing and jet quenching due to parton 
energy loss in dense matter. For definiteness, central collisions are restricted to events 
with impact parameters smaller than 3 fm. 

Hadron spectra at large transverse momentum in pp collisions at collider energy ^/s = 
200 GeV have been measured and are consistent with pQCD parton model calculations 
[0]. The hadron spectra for pp collisions at RHIC energy are expected to be almost the 
same as pp in central rapidity region. Nuclear effects like multiple parton scattering and 
modification of parton distributions in nuclei will modify the hadron spectra in p + A as 
data have shown. These observed nuclear effects at energies -y/i =20 - 40 GeV are well 
accounted for by including the broadening of initial pt due to multiple parton scattering 
and the measured nuclear modification of quark distributions in the pQCD parton model 
[0]. Extrapolation to RHIC is straightforward with a maximum of uncertainty of about 
25% for Pt spectra (normalized to the target nucleons A) at around Pt=4 GeV/c. The 
uncertainty mainly comes from nuclear shadowing of gluon distribution in nuclei which 
only becomes important at or above RHIC energy. Such uncertainties are significantly 
larger in Au + Au collisions as shown by the upper set of curves in Fig. ^ where the ratios 

D . X daAA/dyd^PT 

Raa{Pt) = jr; ^ ... , , „ (5) 

are plotted and (A^binary) is the averaged number of binary A^A^ collisions in A+A collisions 
from the nuclear geometry. Two different parametrizationsQ,^] of the gluon shadowing 
are used in the calculation. If the leading partons suffer medium-induced energy loss, 
the spectra ratios will be suppressed as shown by the lower set of curves in Fig. |[ So 
far, A + A collisions at SPS energy have not shown any indication of energy loss it 
is difficult to predict quantitatively the parton energy loss dE/dx and that is where the 
biggest uncertainty in pf spectra comes from. However, any suppression of the spectra 
at large px will be a clear indication of parton energy loss which will be a direct evidence 
of early parton thermalization. 

The estimate of dNch/drj in A + A collisions at the RHIC energy is much more uncertain 
than the high px spectra because of the dominance of soft physics involved. Here we use 
HIJING model which incorporate both string-like soft production mechanism and 
minijet production. We include both the effect of parton shadowing and jet quenching 
(parton energy loss or partial thermalization of hard partons). The formation of QGP 
or the early thermalization of soft partons and the subsequent expansion will add more 
uncertainty to the estimate of the final observed dNch/drj. 

Shown in Fig. ^ are estimates of dNct/dr] from HIJING with different parametrizations 
of shadowing. Uncertainties due to effect of jet quenching are also shown with or without 
parton energy loss. Including all these effects, there are an uncertainty of about a factor 
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of 2 in the final dNch/drj. Taken the value of dNch/drj = 2.5 in pp collisions at the RHIC 
energy and assuming the A-scaling at RHIC energy is the same as at SPS {PbPb/pp ~ 
250), one has the lowest limit of dNch/dt] = 550 consistent with the HIJING's lowest 
estimate. The clarification of jet quenching from high pj- spectra and of the shadowing 
from p + A data will help to narrow down the uncertainty and help us to understand the 
experimental value of dNch/drj. 




Figure 4. Ratio of hadron spectra in central 
Au + Au and pp collisions as defined in Eq. |^. 



Figure 5. dNch/dr] in central Au + 
Au collisions. EKS § and HIJ || 
parametrizations of parton shadowing are 
used. dE/dx = 1 GeV/fm is used for jet 
quenching scenario. 
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1.4. S.A. Bass: VNI+UrQMD 

Transport theory offers the unique possibihty to cast the entire time evolution of an 
ultra-relativistic heavy-ion reaction at RHIC energies - from its initial state to freeze-out 
- into one consistent framework. In microscopic transport models the full space-time 
evolution of all microscopic degrees of freedom - either all hadrons present in the system 
or alternatively (at higher beam energies) partons - is calculated from the initial state to 
the final freeze-out. At RHIC, the description of the full reaction dynamics from initial 
state to freeze-out needs to incorporate both, partonic and hadronic, degrees of freedom 
explicitly in a consistent fashion. 

In this contribution fl]], a combined microscopic partonic/hadronic transport scenario is 
introduced: the initial high density partonic phase of the heavy-ion reaction is calculated 
in the framework of a parton cascade (VNI) 0, using cross sections obtained in the 
framework of perturbative QCD. The partonic state is then hadronized via a configuration 
space coalescence and cluster hadronization model and used as initial condition for a 
hadronic transport calculation using the Ultra- Relativistic-Quantum- Molecular-Dynamics 
(UrQMD) 0] hadronic transport approach, which has been extensively tested in the AGS 
and CERN/SPS energy regimes. 

The upper left frame of Fig. ^ shows the time evolution (in cm. time, tc.m.) of the rapid- 
ity density dN/dy of partons (i.e. quarks and gluons) and on-shell hadron multiplicities 
at \yc.m.\ < 0.5. Hadronic and partonic phases evolve in parallel and both, parton-parton 
as well as hadron-hadron interactions occur in the same space-time volume. The overlap 
between the partonic and hadronic stages of the reaction stretches from tc.m. ~ 1 fm/c 
up to tc.m. ~ 4 fm/c for the mid-rapidity region. The calculation indicates that this 
overlap occurs not only in time but also in coordinate space - partonic and hadronic 
degrees of freedom occupy the same space-time volume during this reaction phase [^]. 
Hadronic resonances like the A(1232) and the p(770) are formed and remain populated 
up to tc.m. ~ 15 — 20 fm/c, indicating a considerable amount of hadronic rescattering. 

Rates for hadron-hadron collisions per unit rapidity at Ucm. are shown in the left lower 
frame of Fig. ^. Meson-meson and meson-baryon interactions dominate the dynamics of 
the hadronic phase. Due to their larger cross sections baryon-antibaryon collisions occur 
more frequently than baryon-baryon interactions. 

These hadronic rescatterings influence strongly the mrspectra of final state hadrons, 
e.g. negatively charged hadrons and protons, as can be seen in the right upper and 
lower frames of Fig. p. Here calculations with (full circles) and without (open squares) 
hadronic rescattering are shown: both, the low-m^ as well as the high-mt domains of 
the proton spectrum show a depletion due to hadronic rescattering. The proton yield is 
depleted by baryon-antibaryon annihilations which manifests itself mostly in the low-mt 
domain. For h~ hadronic interactions lead to a depletion in the high-mj area of the 
spectrum and to an enhancement in the low-m^ domain. This effect may significantly 
affect signatures of partonic in-medium physics, like jet-quenching. However, the amount 
of hadronic rescattering crucially depends on the treatment of hadronization: if hadrons 
require an additional formation time tj ~ '-y/Mhad after having coalesced from colored 
constituents into color- neutral hadrons before being able to interact, then rescattering 
effects are strongly suppressed (full triangles). 
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VNI + UrQMD 
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Figure 6. Top left: time evolution of the parton and on-shell hadron rapidity densities 
at cm. for central (b < 1 fm) Au+Au collisions at RHIC. Bottom left: Rates for hadron- 
hadron collisions per rapidity at cm. . Right: transverse mass spectra for negatively 
charged hadrons (top) and protons (bottom) with and without hadronic rescattering. 
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1.5. A. Dumitru: Hydrodynamics at RHIC 

I summarize results for Au+Au at RHIC obtained witliin hydrodynamics (also combined 
with a microscopic kinetic model for the late stage). Refs. provide more figures and 
papers, where also references to previous work by others can be found. 

The hydrodynamical description of heavy-ion collisions is very appealing because for a 
given equation of state the (drastic) assumption of local equilibrium uniquely determines 
the solution within the forward light-cone from the specified initial conditions. Here, an 
EoS with first-order phase transition between a hadronic phase and a quark-gluon phase 
is employed, cf. [|| for details. A net baryon and entropy rapidity density of dNB/dy = 25, 
dS/dy = 5000, are employed, as might be appropriate for Au+Au at RHIC. The spectra 




2 3 4 1 2 3 

p,. (GeV) p., (GeV) 



Figure 7. p^-spectra of various hadrons at hadronization (left) and at T = 130 MeV 
(right) from hydrodynamics w/o final resonance decays (Au-|-Au at RHIC, y = 0, b = 0). 

of the hadrons on a given hypersurface can be obtained from the well-known Cooper- 
Frye formula. The p-p-spectra at hadronization and on the T = 130 MeV isotherm are 
shown in Fig. |^. One observes the increasing "stiffness" with the particle mass, i.e. 
(Pt) increases. Also, due to the "hadronic explosion" after hadronization collective flow 
increases substantially if ideal flow persists even deep in the hadronic phase (which it 
doesn't, see below). The work performed by the isentropic expansion (mostly by that of 
the QGP) reduces the transverse energy per unit of rapidity, dEx/dy, by 40% (initially it 
is 1.2 TeV). 

In later stages of the reaction the fluid becomes rather dilute and ideal flow breaks 
down (of course, there may be no region of validity at all). For example, dissipation 
will increase together with the mean-free paths. Furthermore, one can not determine the 
freeze-out hypersurface(s) of the various hadron species from ideal hydrodynamics. 

A solution to these problems is to chose a suitable hypersurface were to switch from fluid 
dynamics to a more detailed kinetic theory (e.g. the Boltzmann equation) which explicitly 
accounts for the various processes between the particles of the fluid. One can then calcu- 
late self-consistently the freeze-out domains of the various elementary processes as deter- 
mined by the local expansion rate and the cross-sections of those interactions, and predict 
hadron momentum spectra, two-particle correlators, etc. Since the equilibrium-limit of 
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the microscopic transport theory is just hydrodynamics, the switch can be performed on 
any hypersurface within the region of vahdity of fluid dynamics (the energy-momentum 
tensors and conserved currents match on the hypersurface). 

Calculations have been performed where the switch was done on the hadronization 
hypersurface. The average collective flow velocity (of it, K, N, A) increases from 0.35 
at hadronization to 0.4-0.45 at freeze-out. The late hadronic "explosion" is found to 
be weaker than at SPS, i.e. kinetic freeze-out at RHIC is closer to hadronization and 
the system cools down less far. In particular, S and fl baryons, which have small cross 
sections, decouple shortly after being hadronized and do not gain substantial flow during 
the hadronic stage. Therefore, the inverse slopes of their m-r-spectra are predicted to be 
smaller than on the T = 130 MeV isotherm, cf. Fig. M. Another interesting observation is 
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Figure 8. Inverse slopes of various hadron species (Au+Au at RHIC, y = 0, b = 0) for 
pure hydrodynamics (on the T = 130 MeV isotherm) and for the Hydro-|-UrQMD model. 



that chemical reactions among the most abundant hadron species (tt, K, N) are difficult 
after hadronization, the respective dN/dy^s change by < 20% (besides the contribution 
from resonance decays). This is because the large expansion rate suppresses inelastic 
processes, which typically have larger relaxation times than elastic scatterings. 
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1.6. H.J. Drescher: NEXUS at RHIC 

neXus is a model to describe high energy interactions from e^e~ annihilation up to 
nucleus-nucleus collisions. The basic physics features (like fragmentation, parton cascad- 
ing, etc.) are tested and parameters are fixed considering the simplest reaction possible. 
This reduces possible sources of errors and uncertainties. 

Basic interaction mechanism in neXus are parton ladders, calculated as a hard scat- 
tering with corresponding initial and final state evolution of partons. The initial state 
radiation is based on a forward evolution algorithm in order to treat consistently energy 
sharing in multiple scattering. 
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Figure 9. Symbolic representation of the total cross section in NEXus. 



It is important to notice that multiple collisions happen simultaneously and momentum 
has to be shared among all collisions according to the different contributions to the total 
cross section. This leads to a complicated formula with a sum over all possible number 
of cut and uncut parton ladders. To choose one of these configurations in a Monte-Carlo 
framework is a highly dimensional problem. This is achieved with techniques coming 
from statistical physics, the Metropolis algorithm and Markov chains. At this stage 
the algorithm uses some approximations, but a future version of neXus will treat this 
properly. 

After primary interactions and particle productions via string decay the rescattering 
of particles is done. For this, particles follow classical trajectories in a Lorenz invariant 
frame. According to cross-sections, two particles do an elastic or inelastic scattering. If 
there are more than two particles close to each other a quark droplet is formed which 
evolves according to a Bjorken scenario. Dropping below some critical energy density, it 
decays according to n-dimensional phase space decay. This reproduces well experimentally 
observed behavior like strangeness enhancement but is based on a theoretically less firm 
ground. 

neXus calculations for gold on gold at 200 GeV center of mass energy and an impact 
parameter less than 2 fm show a charged particle multiplicity of about 1100 particles per 
rapidity unit in the mid-rapidity region (Fig. [To|(a)). 
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Figure 10. (a) Multiplicity of charged particles (full line) compared to SPS lead- lead 
(dashed line), (b) Transverse energy dEt/dy of all particles (full line) compared to SPS 
(dashed), (c) Rapidity of net protons (full line) compared to SPS results (dashed line), 
(d) (full line) and K~ (dashed) and the corresponding SPS results (light dashed lines). 



The transverse energy dEt/dy\y=Q of hadrons is about 850 GeV which leads to an energy 
density of about 7 GeV/fm^ (Figure p!0|(b)). These two observables should have an error 
bar of about 20% due to approximations in the primary interactions and uncertainties in 
the rescattering formalism. 

More significant should be the result for the rapidity distribution of net protons since 
it is quite insensitive to secondary interactions. We find a great transparency with a 
multiplicity of about 5 per rapidity unit at y = 0, which means very little stopping power 
compared to the SPS results (Fig. |rU|(c)). 

Kaons are produced with a multiphcity of about 60 per rapidity unit in the central 
region (Fig. |lO|(d)). 

More plots can be found in 0. 
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1.7. H. Sorge: RQMD at RHIC (see in this proceedings) 
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1.8. M. Bleicher: URQMD at RHIC 

As a tool for our investigation of heavy ion reactions at RHIC the Ultra-relativistic 
Quantum Molecular Dynamics model (UrQMD) is scrutinized. UrQMD is a microscopic 
transport approach based on the covariant propagation of constituent quarks and diquarks 
accompanied by mesonic and baryonic degrees of freedom. It includes rescattering of 
new and produced particles, the excitation and fragmentation of color strings and the 
formation and decay of hadronic resonances. For further details about the UrQMD model, 
the reader is referred to Ref. 0. Note, that these proceedings and the used UrQMD 
version is based on 0. 

The UrQMD model has been applied very successfully to explore heavy ion reaction 
from below AGS energies (10 AGeV) up to the full CERN-SPS energy (160 AGeV). This 
includes detailed studies of particle abundances, -spectra, photonic and leptonic probes, 
J/\E''s and event-by-event fluctuations [Q. 

Let us now start to tackle the relevant questions given to us by the upcoming of the 
RHIC directly: 

• Is a model like UrQMD applicable to AA reactions at RHIC energies? 

• What is the predicted amount of baryonic stopping achieved at RHIC? 

• How many particles will be produced? 

• How strong will the transverse expansion be? 

Let me start by answering the first question: It has often been claimed that dual string 
models fail to describe data above a certain center of mass energy. Let us take a look if this 
region is reached at RHIC: Comparing the UrQMD predictions to rapidity distributions 
in He-|-He at a/s = 31GeV (ISR) leads to a very good description, if we increase to 
energy further UrQMD starts to deviate from the measured multiplicity data by 35% at 
^/s = 200GeV. 

However the collision spectra of individual collisions shows that less than 20% of all 
reactions belong to this high energy regime {^/s > lOOGeV). Thus, the remaining 80% 
of reactions are 'low energetic' and therefore well treatable by string dynamics. 

It has been claimed recently, that exotic mechanisms (Baryon junctions) need to be 
invoked to understand the baryon number transport at SPS and RHIC. In contrast to 
these approaches the UrQMD model mainly applies quark model cross sections to the 
subsequent scattering of constituent (di-)quarks in combination with a small diquark 
breaking component 

The stopping power obtained with this approach is rather strong. The net-proton 
distribution is shifted by more than one united in rapidity, resulting in approximately 12 
net-protons at midrapidity. The stopping of hyperons is even stronger (2 units in rapidity, 
net-A -|- E) leading to a large hyperon density in the central region. 

The particle production predicted by the UrQMD approach reaches more than 1100 
Pions (750 of them are charged Pions) at midrapidity with approx. 100 additional charged 
Kaons at midrapidity. 

This amount of charged particle abundances at midrapidity is on the lower bound of the 
expected multiphcity at y = which reaches from 600 to 2000. It is interesting to note 
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that the particle yield is similar to the presented RQMD results and strongly deviates 
from the reported results obtained by pQCD dominated models. 

Let us turn to the transverse expansion created in the UrQMD model. Since the early 
UrQMD dynamics is based on string degrees of freedom, the newly created quarks are not 
allowed to interact until they have undergone their coalescence into hadrons (typically this 
needs Ifm/c in the local rest frame of the coalescing quarks). Due to the large Lorentz 
7-factor, this leads to a vanishing pressure in the initial reaction zone. This behaviour 
is clearly visible in the distribution of the mean transverse momentum as a function of 
rapidity. A plateau in the {p±){y) is observed for the newly created mesons in the central 
rapidity region - the transverse momenta aX y = are similar to their pp values. In the 
case of the Protons a clear dip in the the {p±){y) at midrapidity is predicted. 

Let me conclude that the semi-hadronic UrQMD model has been applied to Au+Au 
reactions at RHIC. The use of such an effective model for RHIC may be disputable. 
However, rigorous QCD predictions (i.e. results which can be directly compared to data 
to falsify or support this approach) are not available or are even contradicting each other. 
Therefore it is of utmost importance to use a reliable phenomenological model, like the 
UrQMD, which has successfully described a large body of measured data at AGS and 
SPS to investigate the dynamics encountered at RHIC. 
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1.9. B. Zhang: HIJING+ZPC+ART 

To study heavy ion coUisions at the Relativistic Heavy Ion Colhder (RHIC), we have 
developed a muhi-phase transport model that includes both initial partonic and final 
hadronic interactions. Specifically, the parton cascade model ZPC [|| is extended to 
include the quark-gluon to hadronic matter transition and also final-state hadronic inter- 
actions based on the ART model [Q. 

Currently, the ZPC model includes only gluon-gluon elastic scatterings with its cross 
section taken to be the leading divergent cross section regulated by a medium generated 
screening mass that is related to the phase space density of produced partons. The 
input of the parton phase space distribution to the ZPC model is obtained from the 
HIJING model 0, which includes both hard scatterings via the PYTHIA routines with 
the nuclear shadowing effect and soft interactions using the Lund soft momentum transfer 
model. Partons are produced from these scatterings with a formation time determined 
according to a Lorentzian distribution with a half width given by the ratio of its energy 
to the square of its transverse mass. The positions of formed partons are then calculated 
by using straight line propagation from their parent nucleon positions. 

Once partons stop interacting, they are converted into hadrons using the HIJING frag- 
mentation scheme after a proper time of approximate 1.2 fm. We consider both the default 
fragmentation scheme of treating a diquark as a single entity and a modified one which 
allows the formation of diquark-antidiquark pairs, that fragment into both BMB and BB 
with probabilities of 80% and 20%, respectively. 

For the hadron evolution, we have used the ART model that includes both elastic and 
inelastic hadron scatterings. Multiparticle production is modeled through the formation 
of resonances. Since the inverse double resonance channels have smaller cross sections 
than those calculated directly from the detailed balance relation, we have adjusted the 
double resonance absorption cross sections to improve the model. 
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Figure 11. Transverse energy rapidity distribution for central (b=0) Au+Au collisions 
at RHIC. 
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Hadron transverse energy rapidity distributions for RHIC Au+Au central (b=0) col- 
lisions are shown in Fig. |1T[ We see that parton evolution and fragmentation lower the 
distribution by about 100 GeV per unit rapidity while hadron evolution further decreases 
it by around 50 GeV per unit rapidity. This shows an observable effect of final state 



interactions. Fig. gives the baryon and meson rapidity distributions. Net baryon 
distribution has a large rapidity shift due to the modified fragmentation scheme. Many 
antiprotons are seen to survive the absorption in the hadronic matter, leading to a value 
of about 10 at central rapidities. The final meson central rapidity distribution shows a 
distinctive plateau structure. Results using the default HIJING show a similar distribu- 
tion except that the central rapidity density is higher. Also shown in the figure is the 
distribution of kaons produced from both string fragmentation and hadronic production. 
The latter gives significant enhancement of the kaon production. 
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1.10. W. Gassing: VNI+HSD 

Intuitively one expects that the initial nonequilibrium phase of a nucleus-nucleus colli- 
sion at RHIC energies should be described by parton degrees of freedom whereas hadrons 
are only formed (by 'condensation') at a later stage of the reaction which might be a 
couple of fm/c from the initial contact of the heavy ions. Thus parton cascade calcula- 
tions - including transitions rates from perturbative QCD - should be adequate for 
all initial reactions involving a large 4-momentum transfer between the constituents since 
QCD is well tested in its short distance properties. The question, however, remains to 
which extent the parton calculations can be extrapolated to low , where hadronic scales 
become important, and under which conditions (temperature, quark chemical potential) 
the expanding system of strongly interacting partons is more appropriately described by 
an expanding, but interacting hadron gas. 

The practical question is, however, if nonequilibrium partonic and hadron/string models 
can be distinguished at all, i.e. do they lead to different predictions for experimental 
observables? We will thus compare the predictions from the parton cascade VNI with 
those from the HSD transport approach , which involves quarks, diquarks, antiquarks, 
antidiquarks, strings and hadrons as degrees of freedom. 

We start with collisions at a/s = 200 GeV. The calculated results for the proton, 
TT"*" and rapidity distributions in the cms are shown in Fig. |l^ (upper part) for both 
models, which are denoted individually by the labels VNI (dashed histograms) and HSD 
(solid histograms). On the level of collisions we find no considerable differences between 
the two kinetic models with respect to the rapidity distributions for p, vr's and K. 
This also holds for the pt spectra. 

We directly step towards central collisions (6 < 2 fm) for Au + Au at ^/s = 200 GeV. The 
calculated results for the proton (here p—p), antiproton, vr"*" and rapidity distributions 
in the cms are shown in Fig. ^ (middle part) for both models, which are denoted again by 
the labels VNI (dashed histograms) and HSD (solid histograms). Here the hadron-string 
model shows a larger stopping than the parton cascade (l.h.s.) and a flatter distribution in 
rapidity of antiprotons than the partonic cascade. The pion and kaon multiplicities turn 
out to be roughly the same, but the results from the parton cascade are more strongly 
peaked around midrapidity as those from the HSD approach. In physical terms the larger 
stopping of the HSD approach and the additional production of pions and kaons (relative 
to pp collisions) stems from secondary and ternary reaction channels in the hadronic 
rescattering phase, which are quite abundant since the meson densities achieved are very 
high. The narrow width of the antiproton and meson rapidity distribution (as compared 
to pp collisions) for VNI is due to an approximate thermalization of the partonic degrees 
of freedom; in the hadron-string scenario essentially 'comover' scattering occurs with a 
low change of the meson rapidity distribution. Thus the meson rapidity distributions are 
practically the same as for pp collisions. Also note that at midrapidity the net baryon 
density ^ Np — Np is practically zero, however, even at midrapidity at lot of baryons 
appear that are produced together with antibaryons. Thus also mesons (especially cc 
pairs) will encounter a lot of baryons and antibaryons on their way to the continuum. 

In order to demonstrate that the higher baryon stopping, kaon enhancement and widen- 
ing of the pion rapidity distribution is essentially a rescattering effect, we have extended 
the parton cascade VNI by the hadronic rescattering phase as described by HSD 0. The 
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Figure 13. The rapidity spectra of hadrons for pp (upper part) and central collisions of 
Au + Au at y/^ = 200 GeV within VNI, HSD and VNI + HSD (see text). 



results of these calculations (VNI + HSD) for central Au + Au collisions are displayed 
in the lower part of Fig. |l^. It is worth to note that VNI+HSD almost give the same 
rapidity distributions as HSD which already includes the final state interactions. 

We also have performed calculations for low mass dileptons and high mass dimuons , 
where especially the J/\l/ and ip' peaks are of interest. Without explicit representation we 
note that in the dimuon spectra - in the 'comover' absorption scenario - the peak 
vanishes almost completely whereas the J/\l/ peak is suppressed by ~ 90% in central Au 
+ Au collisions due to the high meson densities achieved at these energies. 
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2. HADRONIC FLAVOR OBSERVABLES 

2.1. J. Zimanyi: Quark Coalescence 

There is a possibility, that at RHIC a nice quark gluon plasma will be produced in 
the collision. This plasma will expand and cool in the next phase. Thus it may happen, 
that just before the hadronization the same type of matter will be created as in SPS 
experiments. It is assumed to consist of constituent quark-antiquark matter (CQM) 0,0]. 

Then the hadronization may proceed via quark coalescence process as described by the 
ALCOR model just as in the case of SPS. Thus it is worthwhile to make predictions for 
RHIC assuming coalescence mechanism. 

The nonlinear ALCOR coalescence model was created for situations, where the subpro- 
cesses are not independent, they compete with each other. In this model the coalescence 
equations relating the number of a given type of hadron to the product of the numbers 
of different quarks from which the hadron consists reads as: 

NA = CA-b,-b,-b,-Ng-N,-Ns N^=Cj^-k!-k!-h-N^-N^-Ns 
Ns = Cs-bg-bs-bs-N,-N,-Ns = ■ b^ ■ h ■ h ■ ■ ■ Nj 

Nn = Cn-bs-b,-bs-N,-N,-N, = . . . . . Nj ■ (6) 

= C^.bg-bg- Ng- Ng 

Nk = CK-bg-bj- Ng- Nj 
Nk = C-i^-h-bs-Ng-Ns 

= Crj-bs-bs-Ns-Ns (7) 

Here the normalization coefficients, bi, are determined uniquely by the requirement, that 
the number of the constituent quarks do not change during the hadronization — which 
is the basic assumption for all quark counting methods: 

Ns = 3-Nn + 2-Ns + l-NA + l-N-i^+l-N^ 
Ng = 3- Np + 1- Ns + 2- Na + 1- Nk + 1- N^ 

Ng = 3- Np + 1- + Nj+1- . (8) 

In eq. N^^ is the number of directly produced pions. (Most of the observed pions are 
created in the decay of resonances.) Substituting eqs. - 0) into eq. (|]) one obtains 
equations for the for normalization constants. These constants are then given in terms 
of quark numbers and Ci factors. However, one can predict some relations even without 
solving the set of nonlinear equations, eq. (^ 0,1^ • Thus one obtains the interesting 
relations for the antiparticle - particle ratios: 

A p E A n E 

A p ^ A \l ^ 

where D = K~^/K~ (see details in Ref. f^). At SPS energy the measured value is 
D = 1.8 ± 0.2 P] which agrees within the experimental errors with D = 1.95 obtained in 
ALCOR. Our prediction for RHIC energy is D = 1.2. 
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Figure 14. Particle ratios in Pb + Pb collisions calculated for SPS (open circle) and 
predicted for RHIC (open square). Experimental data (filled circle) are from Refs. 
In Fig. |l^(a) the ratios for SPS and RHIC are practically the same. 



Some particle ratios are shown for the SPS and RHIC energies in Fig. 14. For the 
RHIC calculations the same parameters were used as for the SPS, except the number of 
produced quark antiquark pairs {Nqg = 834, A^^^ = 92 at SPS; Ngg = 4200, A^^^ = 462 for 
RHIC). This number was determined by the requirement, that the same h~ multiplicity 
should be obtained as from HIJING. 
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2.2. J. Stachel: Fireball Model Predictions 

Predictions for hadron yields in the central region at RHIC can be based on a 
thermal model describing successfully the data obtained in present fixed target experi- 
ments at the AGS 0] and SPS [^,|| for central collisions of Si/S and Au/Pb beams with 
heavy targets (Au,Pb). The statistical model treats the system as a grand canonical en- 
semble and has two free parameters, a temperature T and a baryon chemical potential 
IJ,B- Interactions of the particles are taken into account in form of an excluded volume 
correction corresponding to repulsion setting in for all hadrons at a radius of 0.3 fm 
The hadron yield ratios resulting from this model agree particularly well with the most 
complete set of data obtained sofar, the one for central Pb + Pb collisions at the SPS. 
Putting these results in perspective vis-a-vis the expected phase boundary between the 
quark-gluon-plasma and the hadron gas the hadrochemical freeze-out points are where one 
expects the phase boundary and it it suggestive that hadron yields are frozen at the point 
when hadronization is complete PJ^. This suggests for RHIC a hadrochemical freeze-out 
temperature of about 170 MeV, the same as found at the SPS. The chemical potential 
is going to be small; here a value of 10 MeV is used as an upper limit. Strangeness and 
charge conservation require then values of fig and /x/g of 2.5 and -0.2 MeV. 



The yield ratios as predicted for RHIC using these parameters are shown in Figure 15 



The predictions are shown together with a best fit from the same model to data for 
central Pb + Pb collisions at the CERN SPS. One notices drastic differences as compared 
to present data, in particular in the yield ratios of antibaryons to baryons which are 
predicted at RHIC to be about 1 (see Fig. 115]) due to the small chemical potential. 
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Figure 15. Hadron yield ratios from a thermal model (see text) 

In order to predict absolute yields one has to estimate the volume per unit rapidity at 
the time when hadronization is complete. Starting from an initial temperature of Tj = 
500 MeV at a time = 0.2 fm/c and using a transverse expansion in the plasma of jStp = 
0.16 this volume is 3600 fm'^ at a freeze-out temperature of 170 MeV. A selection of the 
predicted hadron yields are shown in Table |l]. The model was used to compute yields for 
193 different hadron species in total. The typical systematic error of this prediction, based 
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on the freeze-out parameters, is judged as follows. The chemical potential is bounded at 
a lower limit of 0. Yields for /if, of 1 MeV were found to not be noticeably different. One 
finds small differences in the absolute yields of baryons with respect to antibaryons. The 
corresponding results are given in column 3 of Tab. |I] where they differ from the values for 
/is = 10 MeV. Further, present SPS data put a lower limit of 160 MeV on the freeze-out 
temperature. Typical results for this lower temperature are given in column 4 of Tab. |I]. 
Generally particle yields are reduced by about 10 % (as indicated for pions in Tab. |I]), for 
the heaviest particles the reduction is 35 % {fl, d). 

Table 1 

Hadron yields per unit of rapidity near mid-rapidity 
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2.3. J. Rafelski: Strange Hadrons from QGP at RHIC 

(Strange) hadronic particle abundances and spectra have been obtained assuming 
a thermally equilibrated, chemically non-equilibrated deconfined QGP source undergoing 
a sudden freeze-out at RHIC. The theoretical framework of this approach is presented in 
some detail in the contribution by Jean Letessier ||2|. 

From the study of the SPS data we deduce the universality of physical properties of 
hadron chemical freeze-out , occurring practically at the same condition as the kinetic 
freeze-out. The differences in system dependent particle m^-slopes are understood to be 
a consequence of differences in collective flow in a deconfined QGP source. We believe 
that this remains the general situation at the 10 times higher collision energies reached 
at RHIC. We expect that the QGP break-up temperature Tf^^ ~ 145 ± 5 MeV will 
see a minor upward change to say Tj'™'^ ~ 150 ± 5 MeV. Thus all we need to do in 
order to characterize particle production at RHIC is to adapt the picture of sudden QGP 
hadronization we have developed to the greater longitudinal flow, lower baryon density 
environment we expect at RHIC. 

In our understanding of the freeze-out, there is complete absence of a pure hadronic 
gas, or even a mixed phase in the dense matter evolution. The deconfined QGP state 
evaporates over a few fm/ c in time, during which time it is kept by a balance of dynamical 
processes near to the freeze-out temperature: the evaporation and work done against the 
vacuum cool the surface down, which is internally heated by the explosive energy flows. 
Empirical evidence at SPS suggests that the transverse flow velocity imparted on emitted 
hadrons does not surpass v± ~ c/a/3 = f sound- The other QGP properties imparted on 
hadronic particles within this sudden freeze-out model are: 

Xs = I ^ (s - s) = ; 

7g ~ 1.5-1.8 -v^ entropy enhancement in plasma and gluons converted to quarks; 

7s > 1 -v^ strangeness equilibrates with T > Tj > Tf] 

unknown : Xg baryon density; baryon:energy stopping relation. 

However, A, must be smaller than seen at SPS in S-S (~ 1.4) so the range Xg ~ 1.2 ±0.15 
which we base our predictions on, is nearly certainly correct. 

We thus expect to find, at RHIC, a kinetically equilibrated, but chemically evolving, 
QGP source, with direct hadron emission from a flowing surface, just as we found in the 
case of SPS energy collisions. We expect, based and compared to the Pb-Pb 158^4 GeV 
results: 

• Shape identity of RHIC m± and y spectra of antibaryons p , A , S , since in our 
approach there is no difference in their production mechanism; 

• The m^-slopes of these antibaryons should be very similar to the result we have from 
Pb-Pb 158A GeV since only a slight increase in the freeze-out temperature occurs, and 
no increase in collective transverse flow is expected; 

• Major changes compared to SPS must be expected in rapidity spectra of mesons, 
baryons and antibaryons. 

Work on an extension of our sudden hadronization model to include the expected sub- 
stantial longitudinal flow at RHIC is underway. We report here preliminary results that 
deserve immediate attention. Strangeness is here an extremely interesting observable for: 
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Table 2 

Sample of hadron ratios expected at RHIC, see text explanation of the model parameters 



Ratios RHIC-1) RHIC-2) 



H/A 


0.18 




0.02 


0.14 




0.02 


H/A 


0.25 




0.03 


0.19 




0.03 


A/A 


0.49 




0.15 


0.49 




0.15 




0.68 




0.15 


0.68 




0.15 


rj/H 


0.14 




0.03 


0.10 




0.03 


n/H 


0.20 


± 


0.03 


0.15 


± 


0.03 


n/n 




1. 






1. 




(n+o) 

(H+H) 


0.17 


± 


0.01 


0.12 


± 


0.01 


(H+H) 
(A+A) 

K+/K- 


0.21 




0.02 


0.17 




0.02 


1.35 




0.25 


1.34 




0.25 


P/P 


2.9 




1.5 


2.9 




1.5 


A/p 


2.4 




0.3 


1.8 




0.3 


h-/B 


4.6 


-2+10 


4.6 ■ 


-2+10 



RHIC rapidity distribution of all (anti)baryons assuming flat net-baryon rapidity distri- 
bution 



dN/dy 


Baryons 


Antibaryons 


protons 


25 (p) 


8.5 ip) 


A 


42 (A) 


20 (A) 




7.5 (H) 


5 (H) 


n 


1 (n) 


1 (n) 



1. Gluons make strangeness effectively in the early stage of the collision and thus at 
SPS strangeness is more central than baryon number. It is possible that this will be even 
more accentuated at RHIC leading to 'strangeness dominance' at central rapidity. 

2. Central strangeness production will further be enhanced by an increase in initial tem- 
perature compared to SPS. We anticipate at RHIC stunning strange hadron abundance 
anomalies. 

In Table ^ we present two scenarios which differ by the strangeness phase space oc- 
cupancy excess: we take in both cases the chemical freeze-out temperature and ve- 
locity Tf = 150 MeV, v± = 1/a/3 ~ 0.58c and the chemical light quark conditions: 
Xq = 1.2 ± 0.15, 7g = 1.8, as well as QGP strangeness conserving = 1, but, in first 
case, we take 7s = 2 and in second 7^ = 1.5 . To obtain the rapidity densities one of the 
key distributions needs to be known: we use a net baryon density d{B — B)/dy ~ 70, 
which if nearly fiat in Ay ~ ±3, allows to fully account for all participating baryons in 
a zero impact parameter collision. Employing the ratios given as RHIC-1) in Table ^ we 
obtain the results shown in Table |^, which imply hyperon-dominance of the baryon yields. 
RHIC-2) scenario would further enhance the hyperon-dominance. 

The key result we see in this extrapolation of the SPS results to RHIC is that the 
(anti)baryon abundances (rapidity distributions) are strongly strangeness dominated, 
which as we have been arguing in past two decades is the characteristic behavior in 
QGP phase break-up. 
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2.4. S.E. Vance: Baryons, Junctions and Hyperons 

Recently, the baryon junction mechanism was proposed to explain the observed baryon 
stopping at the SPS|l|,|^. The baryon junction appears when writing the QCD gauge 
invariant operator of the baryon, being the vertex which links the three color flux (Wilson) 
lines flowing from the valence quarks. Being a gluonic configuration, the junction can be 
easily transported into the mid-rapidity region in hadronic reactions. There near mid- 
rapidity, as the strings connecting the junction to the valence quarks fragment, a baryon 
is produced being composed of three sea quarks. This gluonic mechanism is also able to 
explain P] a striking, preliminary measured baryon asymmetry observed approximately 8 
units of rapidity away from the proton's fragmentation region in ep collisions at HERA. 

The valence baryon junction mechanism was recently extended and a new mechanism for 
anti-hyperon production was proposed|Q. Like the valence baryon junction mechanism, 
this junction-anti-junction loop {J J) mechanism is also derived from the topological gluon 
structure of the baryon and originates in the context of Regge phenomenology. The {J J) 
mechanism was shown to strongly enhance the anti-hyperon yields in nuclear collisions and 
is needed to provide reasonable anti-hyperon to hyperon ratios. Both of these mechanisms 
were implemented in HIJING/BB, a modified version of the HIJING event generator. 




Figure 16. Predictions for the initial valence proton rapidity distribution (up- 
per two curves) and for the initial valence hyperon rapidity distribution (lower two 
curves) are given for Au-fAu collisions at Ecm = 200 AGeV by HIJING (solid) and 
HIJING/BB (dashed). 
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The effects of tliese meclianisms at RHIC energies were explored using HIJING/BB. 
Estimates for the initial valence p and A rapidity distributions in Au + Au collisions at 
RHIC energies (^s = 200 GeV) for 6 = fm are shown in Figure |16[ HIJING/BB 
predicts approximately twice the initial number of valence protons and five times the 
initial number of valence hyperons of HIJING at mid-rapidity leading to a prediction of 
twice the initial baryon density, p{tq) ^ 2po ~ 0.3/fm^. Estimates for the total p and A 
have also been made (see Figure |17D, where it was found that at these energies, the yields 
of the anti-baryons are only sensitive to the relative string fragmentation probability of 
producing diquarks to quarks, and not to the small cross section for producing J J loops. 
While the J J loops do not significantly effect the absolute yields, they are important in 
producing rapidity correlations between baryons and anti-baryons, such as the A'^~^{uuu) 
and Q~^{sss), which are absent in present baryon pair production schemes. 




Figure 17. Predictions for the initial p and A rapidity distributions are given for HIJING 
(solid lines) and for HIJING/BB (dashed curves) for Au-|-Au collisions at E^m = 200 
AGeV for & = fm. 



In the phenomenology associated with the baryon junction, many new states being 
gauge invariant combinations of junctions, anti-junction, quark and anti-quarks |^ were 
predicted, e.g. a quarkless hybrid glueball (JJ = Mq). The experimental search for 
these states has been inconclusive (presumably due to their large decay widths), with 
only several broad resonances remaining as possible candidates. We propose that the 
effect of large glueball states (e.g. Sq = {J J J J J J)) can be observed at RHIC by seeing 
a non-trivial rapidity correlation between two baryons (e.g. pp rapidity correlations) in 
A^A^ colhsions (see @). 
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3. HBT AND COLLECTIVE FLOW CORRELATIONS 

3.1. S. Pratt: HBT at RHIC 

Prediction of QGP time delay too late for these proceedings. [M.G.: Look for it at 
RHIC in spite of heretical views in sec. 3.3. See D.H. Rischke, M. Gyulassy, NPA 608 
(1996) 479] 

3.2. D. Teaney: Nutcracker Flow and HBT 

Flow in non-central collisions at the SPS has been studied experimentally P] and the- 
oretically and shows sensitivity to the EoS. However, high and low energy collisions 
have different acceleration histories, due to the QCD phase transition. In the mixed phase 
the EoS is soft 0. For AGS/SPS collision energies the matter is produced close to the 
"softest point" and therefore the initial transverse acceleration is small. Using a bag 
model EoS and fixing the initial entropy to the SPS multiplicity, we have simulated PbPb 
collisions at b=8fm with boost invariant hydrodynamics. The matter expands preferen- 
tially in the impact parameter direction (the x direction). Since the matter is produced 
close to the softest point, it retains its initial elliptic shape, accelerating only in the late 
hadronic stages. For higher collision energies, the hot QGP is formed well above the tran- 
sition temperature and the early pressure forces early acceleration. Unfortunately, the 
final radial flow velocity at RHIC is quite similar to the SPS as is the final elliptic flow 
The early quark push has important consequences, however. The resulting velocity 
has a long time, ~ 10/m/c, to change the matter distribution before freeze-out. The 
stiff QGP in the center , with T » , pushes against the soft matter on the exterior, 
with T ?a Tc, producing two shell-like structures which separate in the x direction leaving 
two holes in the y direction. Since the final distorted distribution rather resembles a nut 
and its shell, we call this picturesque configuration the nutshell. In Fig.0(a) we show 
the "nutty" energy distribution in the transverse plane at zero rapidity. No shell like 
structures are seen for an ideal gas EoS. The matter expands elliptically and the final 
distribution resembles a gaussian. 

The early acceleration may be seen with HBT interferometry since the early velocity 
increases the size of the system. From the two pion correlation function, the following 
radii may be extracted by carefully choosing the pion pair momenta 

Rl^ =<x^ > - <x>'^ (9) 
Rly =<y'>-<y>' (10) 

(The average momentum of the pair was taken to be 500 MeV). These radii are shown in 
Fig.|l^(b) for PbPb collisions at b=8 fm, as a function of the pion multiplicity scaled by 
the number of participants to central collisions (not b=8 fm). i?^^. and Ryy are shown for 
a bag model EoS and for a simple resonance gas EoS. For low energies near the left hand 
side of the plot, the two EoS show approximately the same radii, roughly corresponding 
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Figure 18. (a) Typical matter distribution in the transverse plane at mid rapidity, cal- 
culated using boost invariant hydrodynamics with a bag model EoS, for a PbPb collision 
at b=8 fm at RHIC( dNT,/dy = 850 ) . The solid lines show contours of constant energy 
density. The inner and outer dashed lines show temperatures T=140 MeV and T=120 
MeV respectively, (b) HBT radii, R^x and Ryy, as a function of the final pion multiplic- 
ity. For a resonance gas EoS p = .2e, the solid lines and dotted lines show R^x and Ryy 
respectively. For a bag model EoS, the dashed and dashed-dotted lines show Rxx and Ryy 
respectively. 



to the initial elliptic shape of the matter distribution. For a simple resonance gas EoS 
the HBT Radii show little energy dependence, while for an EoS with the phase transition 
the Radii increase steadily with beam energy. The rapid increase of Rxx can be under- 
stood qualitatively. The nutshells separate at higher collision energies and leave a large 
homogeneous region in the center. The constancy of Rxx for an ideal gas EoS reflects 
the energy gradients in the final distribution. For an EoS with a phase transition, the 
source function is "box" like and the "box" increases in size, while for an ideal gas EoS 
the source function is more gaussian. 

Other signatures of the nutcracker flow were presented in and similar qualitative 
features were found by i,|lO|. 

Acknowledgments: Work supported by U.S. DOE, under Contract No. DE-FG02- 
88ER4038. 
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3.3. B.R. Schlei: The QGP Stall or Not 

In 1996, D.H. Rischke and M. Gyulassy [l|] proposed to use the time-delay signal in 
Bose-Einstein correlations (BEG) as a signature for the possible formation of a quark- 
gluon plasma (QGP) in relativistic heavy-ion collisions. In particular, they suggested to 
measure the ratio of the transversal interferometry radii. Rout/ Rside, which can be obtained 
by fitting experimental BEG data with a parametrization introduced by G. Bertsch 
et al. in 1988. The transverse radius parameter Rout has compared to the transverse 
radius parameter Rside an additional temporal dependence which should be sensitive to 
a prolonged lifetime of a fireball, in case a QGP was formed in a relativistic heavy-ion 
collision. This latter phenomenon is known under the term "QGP stall" . In the following, 
I shall explain why I believe that the time-delay signal in Bose-Einstein correlations is 
not a good signature for the possible formation of a QGP. 

Let us consider two quite different equations of state (EOS) of nuclear matter within a 
true relativistic hydrodynamic framework (i.e., HYLANDER-G) 0. The first one, EOS-I, 
has a phase-transition to a QGP at Tc = 200 MeV {ec = 1.35 GeV/fw?) @]. The second 
EOS, EOS-II, is a purely hadronic EOS, which has been extracted from the transport 
code RQMD (cf., Ref. |^) under the assumption of fully achieved thermalization. If one 
assumes for each EOS different initial conditions before the hydrodynamical expansions, 
one can fit simultaneously hadronic single inclusive momentum spectra and BEG, which 
have been measured recently by the GERN/NA44 and GERN/NA49 Gollaborations (cf., 
[^,0), respectively. In particular, for the acceptance of the NA44 experiment a ratio 
Rout/ Rside ~ 1-15 was found while using both EOS, EOS-I and EOS-II. Little difference 
was seen in the BEG of identical pion pairs while considering the two different EOS. 

In the following, we shall assume for central Au+Au collisions at BNL/RHIG beam 
energies a set of fireball initial conditions, IG-I, which are similar to those as described in 
Ref. []^. From these fireball initial conditions, IG-I, single inclusive momentum spectra 
have been calculated while using EOS-I and EOS-II in the hydrodynamic expansions. 
We note, that the rapidity spectra of both calculations differ in width and normalization 
significantly . Fig. shows, that the isotherms of the transversely expanding fluids at 
longitudinal coordinate z = Q also differ significantly. Since there will be only one set of 
measured data, we shall fit the calculation using EOS-II to the single inclusive momentum 
spectra of the calculation using EOS-I. In doing so, we find new initial conditions, IG-II 
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Figure 19. Isotherms in steps of AT = 20MeV. The outer contours equal Tf = lAOMeV. 



0. But now the space-time picture of the evolving fireball at freeze-out is again very 
similar to the one using EOS-I with IC-I. 

If one calculates BEG of identical pion pairs or identical kaon pairs, respectively, one 
finds 1^ while comparing the calculations using EOS-I with IC-I and EOS-II with IC- 
II, respectively, no significant differences in the extracted ratios Rout/Rside regardless of 
the pair kinematics under consideration. In particular, the assumption of the PHENIX 
detector acceptance [^j leads to a ratio Rout/Rside ~ 1-65 for both choices of EOS. 

In summary, the larger ratio Rout/Rside at RHIC beam energies appears to be rather 
a consequence of the expected higher energy deposit in the fireball during the heavy-ion 
collision, but it appears not to be an indicator of the present or absent phase-transition 
to a QGP. Of course, more theoretical analysis is necessary, but there is strong evidence, 
that BEG do not provide a good QGP signature, since we do not understand the initial 
state of a heavy-ion colhsion well enough yet. 
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3.4. J.Y. Ollitrault: Elliptic Flow 
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3.5. D. Molnar: Elliptic Flow via Inelastic Parton Cascade 

Transverse energy and elliptic flow are of recent theoretical interest because they signal 
collectivity and thus may be sensitive to a possible phase transition. Also, they are 
relatively easy to measure and therefore will be one of the flrst results at RHIC 

Until thermal and chemical equilibrium are shown to be achieved in relativistic heavy 
ion collisions, a transport theory approach is absolutely necessary. Two forerunners 
of the present study utilize the ZPC elastic parton cascade p, which has colorless gluons 
with a single elastic 2 — 2 process. However, for the question of chemical equilibration it 
is essential to include inelastic processes, too. 

If we insist on having on-shell partons, the simplest set of particle number changing 
processes consists of a 2 — 3 and a 3 — 2 process. By adding these two new processes to 
ZPC, we developed an on-shell inelastic parton cascade, MFC 0. 

To have a well-deflned theoretical model that is independent of the code, we regard 
MFC as a tool for solving the Boltzmann transport equation 



where / is the Lorentz invariant phase space density, the Lorentz invariant W-s are 
proportional to the corresponding 2 — 2, 2 — > 3 and 3-^2 matrix elements, fi = f{pi,x), 
J = J -J^, and g is the degeneracy of the partons. Furthermore, combinations like 

12 — 34 are for pi +P2—P3 —P4, etc., and the arguments of ^^12^34(^15^2 5^3,^4)5 etc., are 
suppressed. 

For simplicity, we took our gluons to be massless. The 2 — 2 and 2 3 differential 
cross sections were taken to be uniform in phase space with theta function cutoffs 6((pi + 
PjY — fi"^) for each possible pair of ingoing and each possible pair of outgoing particles 
in order to prevent coUinear collisions. The leftover momentum dependence of the cross 
sections was flxed by imposing energy independent total crossections for /x^ = 0, while the 
3 — *• 2 cross section was given by detailed balance. The gluon degeneracy was g = 16 {3 
colors) . 

Fig. ^ shows the transverse energy loss in a one-dimensional Bjorken expansion 
starting from a locally thermal distribution. Keeping the total 2 ^ X cross section 
constant, we studied how the energy loss changes with the ratio of the elastic and inelastic 
part of the cross section. We considered three cases: fully elastic, 50% inelastic, and fully 
inelastic. The initial conditions were the same as in Ref. |l|]: tq = O.lfm/c, Tq = SOOMeV, 
dN/dr] = 400, -5 <r] <5, At = lOOfm^ and we averaged over 20 events. 

We can see that the transverse energy at midrapidity is sensitive only to the total 
2 X cross section. With increasing inelasticity there is somewhat less work done 
because opening up the inelastic channels leads to parton production which compensates 
for the decrease of work. There also shows a tendency that until the cross section is 
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Figure 20. Ef evolution for a 1- dimensional Figure 21. Elliptic flow vs. centrality and 
Bjorken expansion inelasticity 



dominantly inelastic the energy loss is rather insensitive to the degree of inelasticity. 

Fig. ^ shows how the elliptic flow parameter V2 changes with centrality and inelas- 
ticity. We took a similar initial condition as Ref. but with dN/drjd'^xt = 3/fm^, 
To = SOOMeV, R = 7fm, —5 < ?7 < 5. We averaged over 20 events. 

As we go from the fully elastic to the fully inelastic case elliptic flow drops, showing 
that the details of the microscopic dynamics are important. This is also supported by 
another study using the RQMD 3.0 event generator. 

The 50% reduction of V2 is in addition to the original 50% drop relative to ideal hy- 
drodynamics given by dissipative corrections at the elastic level 0. The centrality 
dependence is very similar to that from hydrodynamics. 

A study with more realistic gluon cross sections and diffuse initial geometry is underway. 
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4. JETS AND PENETRATING PROBES AT RHIC 

4.1. R. Rapp: Low Mass Dileptons 

Low-mass dilepton spectra in (ultra-) relativistic heavy-ion collisions are expected to 
provide information on the in-medium properties of the light vector mesons in connection 
with signals from chiral symmetry restoration. At RHIC energies, the major experimental 
challenge is a large (combinatorial) background which limits the sensitivity to well-defined 
resonance structures. In this contribution [|l| we will thus address the question to what 
extent p and uo spectral functions are affected in meson-dominated matter (following 
paragraph) and how this reflects itself in space-time integrated dilepton spectra to be 
measured at RHIC (last paragraph). 

Thermal dilepton radiation from a hot and dense medium is governed by the thermal 
expectation value of the electromagnetic current-current correlator, Hem- Up to invariant 
masses of ~ 1 GeV its hadronic part can be accurately saturated with light vector mesons 
(Vector Dominance Model). The pertinent thermal dilepton production rate can then be 
expressed through the vector meson spectral functions, \mDy, as 

^r-i ^ """-"^ ^' ■ '""-''^ ^' ^ f ^' ■ 

where denotes a thermal Bose distribution function and M is the invariant mass of the 
lepton pair. Medium effects in the p propagator are accounted for through modifications 
in the p — vrvr width (most notably from a Bose enhancement of the two-pion states 
characterizing 'stimulated emission') and by collisions with surrounding thermal vr, K 
and p mesons via s-channel resonance formation (tu, Oi, Ki, fi, etc.), see, e.g., Refs. 
At a temperature T=150 MeV the total thermal broadening of ImDp amounts to about 



80 MeV with no significant shift in mass, cf. left panel of Fig. 22 
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Figure 22. p- (left panel) and uj- (right panel) spectral functions (weighted by the corre- 
sponding VDM coupling constant) at fixed 3-momentum in a thermal meson gas. 
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Along similar lines, the uj meson decay width (which in free space consists of a com- 
bination of pn and direct Bvr decays) has been evaluated including both pion- and rho 
Bose factors as well as the in-medium p spectral function. Also, an in-medium selfenergy 
from Tiuj scattering through the 61 (1235) resonance has been inferred from the 70% 
branching ratio of the free hi —>■ uj-k decay. As a result, the finite temperature uj spectral 
function exhibits a broadening of several times its vacuum width, cf. right panel of Fig. |22|. 

Dilepton spectra from p and uj decays are readily calculated by integrating the thermal 



rate, Eq. (|T2D, over the hadronic space-time volume of a central Au+Au collision at RHIC. 
For that we assumed an adiabatic expansion (including a mixed phase at Tc= 180 MeV) 
with an entropy over baryon density of s/ub = 220 together with a charged particle 
multiplicity of dNch/dy ^ 1100 around midrapidity. The 3-volume expansion has been 
modelled in accordance with hydrodynamic simulations. The final 3-momentum inte- 
grated dilepton invariant mass spectra from p, — > e~^e~ decays are displayed in Fig. 123 



The broad p structure can most likely not be disentangled from the large background ex- 
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Figure 23. Dilepton spectra from a thermal fireball through p, — >■ e^e decays in the 
hadronic phase including a schematic acceptance for the PHENIX detector. 



pected in PHENIX. On the other hand, the thermal broadening of the lo (leading to a net 
~ 40 MeV wide signal) may be 'just right', z.e., narrow enough to be possibly detectable 
above the background and the p contribution, but substantially wider than the free a;, 
which allows for a discrimination from (free) lo decays after freeze-out and thus enables 
important insights into its finite temperature behavior. 
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4.2. I. Vitev: Jet Quenching in Thin Plasmas 

One of the predicted new observables in nuclear collisions at RHIC energies (a/s ~ 
200 AGeV) is jet quenching due to gluon bremsstrahlung. We study the detailed gluon 
probability distribution in the case of few {rig = 1, 2, 3) scattering centers which leads to 
the suppression of the high p± tails of the hadronic spectra . 

The Landau-Pomeranchuk-Migdal (LPM) effect first studied for the case of QED, 
plays an important role in the non-abelian gauge theories as well, leading to nontrivial 
results [Q. Unlike in electrodynamics, where the energy of the final electron and photon 
can be directly measured, the energy of a single parton is not measurable. Therefore, in 
order to test the non-abelian energy loss of jets [Q, the detailed form of the angular distri- 
bution, dNg/dy(Pk±, of induced gluon radiation must be known. The angular integrated 
energy loss, dE/dx, is theoretically interesting but insufficient to address the experimen- 
tally accessible question of how much are moderate ~ 3 — 10 GeV pions expected to 
be suppressed relative to their known spectrum in elementary nucleon-nucleon reactions. 

We use the Gyulassy-Wang (GW) model of locally thermally equilibrated plasma [^]. 
In the case of one scattering center there are three relevant diagrams |@] that give the 
gluon conditional probability distribution: 

= Cn^(H' + R{Bl + p{b,)Cl)-R(H-B^cos{t,ouJo)) 
dyd^k± TT^ ^ ^ ^ 

-Rp{h/2) (#-C'icos(tio^io)-2C'i-5icos(tioa;i))) , (13) 
where we have rearranged the terms following ||^ and introduced the following notation: 



^ ^ ^ {k-qi) _ 



Tj V "iiyj- Tj . 

n — -j-J , — — , -Dl — -H — Ui , 



ujq = kj_/2uj = 1/tf , uJi = {k ~ qi)j_/2uj , loiq = loi — ooq ■ (14) 

The phase factors ujq, uji, ujiq have the formation physics built in. The normalized trans- 
verse density profile p(&i) = T(6i)/T(0) eliminates rescatterings of the emitted gluon at 
large angles and keeps consistency with the geometry of the process. R = Ca/Cr is 1 or 
9/4 for gluon or quark jets respectively. 

Similar formulas cam be derived for the case of two and three scattering centers 0. In 
all cases a couple of simple analytical limits can be obtained. 

lim ^ = ^ , lim ^ = ^ (1 + i?) . (15) 

''^^ dyd'^kx dyd'^k± ''^^ °° dyd'^k± dyd'^k± 

Therefore we can make an ansatz, allowing us to consider multiple scattering extrapola- 
tions: 



dyd^k^ dyd^k± ^ ^ \ uj jj^ ' 



From Fig. 24. one can see that the quantitative agreement between the scaled result for 
one scattering center and the exact and much more complicated = 2 result is within 
5% everywhere in the kinematically allowed domain. 
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Figure 24. Scaling of the gluon prob- 
ability distribution with the number of 
scattering centers. Ejet = 50 GeV 



Figure 25. High p± tail suppression of 
the hadronic distribution for partonic 
event Ec.m. = 20 GeV (JETSET). 



Jet quenching occurs even in free space where there are no final state interactions, 
leading to characteristic double-logarithmic probability distribution In practice, the 
large logarithm implies t hat multiple gluon emission must also be considered. This leads 
to a Sudakov form factor for the jet and a probabilistic Alterelli-Parisi parton shower [0. 

The final state multigluon shower can be calculated most readily by one of the many 
Monte-Carlo event generators that encode empirical parton to hadron fragmentation 
functions and thus allow the detailed study of the effect of parton showering on the final 
hadron distributions. We have used the Lund JETSET string fragmentation routine 
to hadronize a high energy partonic event accompanied by gluon bremsstrahlung. Fig. 25. 
clearly shows the quenching effect of the medium induced and modified radiation. 

Considering the significant role that "hard" physics will play at RHIC energies as 
compared to SPS we expect this effect to be roughly on the order of two. 
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4.3. A.M. Snigirev: Monojet Rate at RHIC 

Hard jet production is considered to be an effective probe for formation of super-dense 
matter - quark- gluon plasma (QGP) in future heavy ion collider experiments at RHIC and 
LHC. High parton pair (dijet) from a single hard scattering is produced at the initial 
stage of the collision process (typically, at ~ 0.01 fm/c). It then propagates through 
the QGP formed due to mini-jet production at larger time scales (~ 0.1 fm/c), and 
interacts strongly with the comoving constituents in the medium. The various aspects 
of hard parton passage through the dense matter are discussed intensively [0-1^. In 



particular, the strong acoplanarity of dijet transverse momentum ||T]-§1, the dijet quenching 
(a suppression of high px jet pairs) p and a monojet-to-dijet ratio enhancement p|] were 
originally proposed as possible signals of dense matter formation in ultrarelativistic ion 
collisions. 

In the simple QCD picture for a single hard parton-parton scattering without initial 
state gluon radiation (i.e. when jets from dijet pair escape from primary hard scattering 
vertex back-to-back in azimuthal plane with equal absolute transverse momentum values, 
Pti = PT2) a monojet is created only if one of the two hard partonic jets loses so much 
energy due to multiple scattering in the dense matter that effectively we can detect only 
one single jet in the final state. The monojet rate is obtained by integrating the dijet rate 
over the transverse momentum pt2 of the second (unobserved) jet with the condition that 
Pt2 be smaller than the threshold value pcut (or the threshold jet energy Et = Pent Q)- Then 
rate of dijets i?'^*-?^* with pti,Pt2 > Pcut and monojets R^°^° with pTi > Pcut {pt2 < Pcut) 
in central AA collisions is calculated as integral over all possible jet transverse momenta 
Pti, Pt2 and longitudinal rapidities yi, y2- 

At first in the framework of the simple model |T^ we demonstrate that monojet-to-dijet 



ratio can be related to mean the acoplanarity measured in the units of the jet threshold 
energy, namely 

TDmono ^ I I \ 

^ ,. . ^ oc . 17 

The results of physics simulation have been obtained in the three scenarios for jet 
quenching due to collisional energy losses of jet partons in mid-rapidity region y = |]TT], 
|T^: (i) no jet quenching, (ii) jet quenching in a perfect longitudinally expanding QGP 
(the average collisional energy losses of a hard gluon < AEg >~ 10 GeV, < AEg >= 
4/9- < AEg >), {iii) jet quenching in a maximally viscous quark-gluon fluid, resulting in 
< AEg 20 GeV. Initial state gluon radiation has been taken into account with the 
PYTHIA Monte-Carlo model ^ at c.m.s. energy ^/s = 5.5A TeV. 

Thus we conclude that rescattering of hard partons in medium results in weaker Et- 
dependence of ratio Rjnono j j^dijet ^.^ ^ \Kt\ > /Et- With growth of energy losses the 
ratio we are interested in has a tendency to be constant, what would be interpreted as 
the signal of super-dense matter formation. 
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Figure 26. (Monojet/dijet)/(mean acoplanarity/£^T) ratio as a function of jet energy 
threshold Et- 
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4.4. D.K. Srivastava: Direct Photons and Dileptons at RHIC 

The radiation of single photons from quark matter has recently been estimated to the 
order of two loops This provides a large bremsstrahlung contribution as well as a new 
mechanism of quark-annihilation with scattering. Together these lead to a much large 
radiation of photons than the estimates based on one- loop calculations [|]. These rates 
along with the radiations from the hadronic matter due to the hadronic reactions have 
been used to estimate the yield of photons from Pb + Pb collisions at RHIC with the 
assumption that a thermally and chemically equilibrated quark gluon plasma is formed at 
To = 0.5 fm/c and having Tq = 310 MeV, corresponding to dN.„/dy = 1734 (see Fig. |27|) . 

However, it is rather unlikely that the QGP would be created in a chemically equili- 
brated form. The radiation of photons from an equilibrating plasma |^ due to Comp- 
ton and annihilation processes (alone) and also from the pre-equilibrium stage (i.e., pre- 
thermally equilibrated quark matter) within a parton cascade model |^ are given in 
Fig. H. 



Pb+Pb § SHIC; Hydrodyaamics 
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Figure 27. Radiation of photons from central collision of 
lead nuclei at RHIC energies from the hadronic matter 
(in the mixed phase and the hadronic phase) and the 
quark matter (in the QGP phase and the mixed phase). 
The contribution of the quark matter while using the 
rates obtained by Kapusta et al. and Aurenche et al. 
and those from hard QCD processes are shown sepa- 
rately. 
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Figure 28. Radiation of photons from central colli- 
sion of gold nuclei at RHIC energies from an equili- 
brating plasma where the initial conditions are taken 
from a self-screened parton cascade model. Only Comp- 
ton and annihilation processes are included. The pre- 
equilibrium contribution within a parton cascade model 
is also shown. 



Correlated charm decay presents the largest source of dileptons at RHIC energies, 
though the spectrum of such dileptons will depend sensitively on the extent of thermaliza- 
tion of the charm quarks and D-mesons P] . In case it is absent the upper limit of dilepton 
production can be seen in Fig. E9|, while the radiation of dileptons due to annihilation of 



quarks (alone) from a chemically equilibrating plasma ^ is shown in Fig. ^ 
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Sources of Dileptons; Au+Au®RHIC 




Figure 29. Sources of high mass dileptons. A thermal- 
ized and chemically equilibrated QGP is assumed. The 
correlated charm and bottom decay is estimated assum- 
ing no energy-loss for the heavy quarks and no thermal- 
ization of D mesons. 



Figure 30. Radiation of large mass dileptons from a 
chemically equilibrating quark gluon plasma. 
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4.5. Z. Lin: Open Charm and Drell Yan at RHIC 

One of the signatures for the formation of the quark-gluon plasma in heavy ion collisions 
at RHIC is the thermal dileptons emitted from such a matter. However, dileptons from 
both final heavy meson decays and initial Drell- Yan processes can contribute significantly 
to the continuum background, thus making it difficult to observe the thermal dileptons 
from the quark-gluon plasma. Quantitative studies of dilepton production from these 
processes are therefore essential. 

Although production of charm quarks in proton-proton interactions can be reason- 



42 



ably well described by perturbative QCD, their fragmentation to charm mesons is non- 
perturbative. In order to calculate the final open charm observables, one needs to know 
the intrinsic transverse momentum kf of the partons inside a proton, the K-factor which 
takes into account higher order contributions, and the parameter e when the Peterson 
charm fragmentation function is used. From the fit to E791 data on the angular dis- 
tribution between the charm and anti-charm mesons, kt is found to be about 1 GeV. In 
fitting the charm meson p±_ spectra from E706 and the spectra from E769, Kc — 3 and 
e ~ 0.01 are obtained. For Drell-Yan processes, we simply take the leading-order (LO) 
qq annihilation cross sections with a K-factor, Kdy = 1-7, determined by fitting to E772 
data. Thermal dileptons from the LO qq annihilation from the quark-gluon plasma is cal- 
culated in the parton cascade model ZPC, where the initial gluon and quark distributions 
are taken from the HIJING model. 



In Fig. P^, we show the dilepton invariant-mass spectra for central Au+Au collisions at 
RHIC. The solid, dashed and dot-dashed curves represent dileptons from correlated open 
charm decays, Drell-Yan processes, and LO thermal qq annihilation, respectively. The 
dilepton yield from open charm decays is lower than earlier estimates where the charm 
fragmentation function is taken to he 5{1 — z), which is equivalent to taking e = in the 
Peterson formula. For dileptons with mass below 10 GeV, the open charm contribution is 
higher than that from Drell-Yan processes, which is much larger than the thermal dilepton 
yield from LO qq annihilation from the quark-gluon plasma. 

Our preliminary results thus indicate that to detect the thermal dileptons from the 
quark-gluon plasma does not seem feasible. However, we have underestimated the ther- 
mal dilepton yield. First, only 2 — > 2 partonic processes are included in ZPC, so the 
total number of partons does not increase as the system evolves. Secondly, we have 
not included thermal dilepton production from gluonic processes such as gq — > q-y* and 
99 ~^ I'll* ■ Although cross sections for these two processes are suppressed by as and a^, 
respectively, they could have large effects on the thermal dilepton yield because gluons 
greatly outnumber quarks in the early partonic stage of heavy ion collisions at RHIC. 

Also, results on open charms are only based on information taken from pp collisions. 
In heavy ion collisions, charm mesons are expected to undergo final-state hadronic inter- 
actions. We have recently studied the hadronic scatterings between charm mesons and 
hadrons such as pion, rho, and nucleon and have found that the charm meson spectra 
are significantly affected by these rescatterings. The increase of the invariant mass of 
charm meson pairs due to scatterings in hadronic matter leads to an increase of the in- 
variant mass of dileptons from their decays, and this has been shown to provide a possible 
explanation for the observed large dimuon enhancement in the intermediate-mass-region 
(IMR) of 1.5 — 2.5 GeV in the NA50 experiment at SPS. Such final-state hadronic rescat- 
tering effects are different from those due to initial-state interactions, e.g., an increase of 
the parton intrinsic kt in the nuclei effectively boosts the charm pair in the transverse 
direction, thus causing negligible change of the pair mass spectra and the resulting IMR 
dileptons from their decays. 

At RHIC energies and beyond, partonic interactions of charm quarks in the quark-gluon 
plasma are also possible, such as the energy loss of fast quarks. Even a moderate energy 
loss of 0.5 GeV/fm has been found to suppress by an order of magnitude both the yield 
of high charm mesons and that of large invariant-mass dileptons from their decays. 
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Figure 31. Dilepton spectra from central Au+Au collisions at RHIC. 



For a more reliable evaluation of the dilepton spectra in heavy ion collisions at RHIC, 
one needs to further study the effects due to final-state interactions of both charm quarks 
and charm mesons as well as thermal dilepton production from gluonic processes in the 
quark-gluon plasma. 
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4.6. R. Thews: Mesons at RHIC 

The Be meson is the bound state of be (or be) whose recent detection is the first step 
toward completion of the spectroscopy of heavy quark mesonic states. The b-e states have 
properties that conveniently fill the gap between the J/ip and the T states. Thus it is 
probable that at RHIC the B^ mesons will serve as a probe of deconfined matter We 
find that significant differences arise for B^ formation in deconfined and confined matter 
[0]. Our initial calculations suggest that: 

(a) The rates of normal hadronic production mechanisms at RHIC energies are not 
sufficient to produce a detectable number of Be mesons. 

(b) If a region of deconfined quarks and gluons is formed, the production (and survival) 
rate can be enhanced by several orders of magnitude. 

(c) The observation of B^ mesons at RHIC would signal a source of deconfined charmed 
quarks, and the rate of Be production will be a measure of the initial density and tem- 
perature of that source. 

We note that the study of the b-e sector has the advantage of a long history of potential 
model analysis in the bb and cc sectors. These studies have provided robust predictions 
for the mass and lifetime of the Be states 0, and the recent measurements by CDF are 
consistent with those calculations 0. 

First, we estimate the production rate of Be which one would expect if it results just 
from a superposition of the initial nucleon-nucleon collisions at RHIC. For heavy quark 
production, pQCD calculations for p-p interactions fit present accelerator data and bracket 
the RHIC energy range. Hard Probes Collaboration estimates indicate about 10 cc pairs 
and 0.05 bb pairs per central collision [Q. J/ip and T production involves the use of some 
model, such as the Hard Probes color singlet fits which predict bound state fractions 
of order somewhat less than the one percent level. A similar analysis for Be production 
involves substantially smaller bound state fractions, since the b and c must be produced 
in the same hard interaction, a process of order a^. At RHIC energies, typical values are 
3 — 10 X 10~^ 0, with the uncertainty from the scale choice in the pQCD calculations. 

To get predictions for RHIC, we have looked at two scenarios for the luminosity, a) 
The "first year" case assumes a luminosity of 20 fib~^ with no trigger, b) The "design" 
luminosity assumes 65 Hz event rate with a 10% centrality trigger in Phenix, and uses 10^ 
sec/year. The event rate predictions listed in the Table include both the weak branching 
fraction of the Be plus the dimuon decay fraction for J/ip. Results are also shown for J/ip 
and T production and detection via and the underlying heavy quark production 

which may be useful to make contact with other estimates. One sees easily that in this 
scenario there is no hope of seeing Be's at RHIC. 

Now the principal reason for this work - could deconfinement change the Be production 
rate at RHIC? We have investigated the following scenario: For those events in which a 
bb pair are produced, the small Be formation fraction could be avoided if one utilized the 
10 cc pairs already produced by independent nucleon-nucleon collisions in the same event. 
This can occur if and only if there is a region of deconfinement which allows a spatial 
overlap of the b and c quarks. Also, one would expect additional cc production in the 
deconfined phase during its lifetime, as a result of the approach to chemical equilibration. 
The large binding energy of Be (840 Mev) would favor their early "freezing out" and they 
will tend to survive as the temperature drops to the phase transition value. The same 
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effect for the B mesons and indeed for the Bg will not be so competitive, since these 
states are not bound at the initial high temperatures (or equivalently they are ionized at 
a relatively high rate by thermal gluons). 

To do a quantitative estimate, we first utilized the quarkonium break-up cross section 
based on the operator product expansion [0 , to calculate the dissociation rate of bound 
states due to thermal gluon collisions. We then used detailed balance for the formation 
rates, and calculated the equilibrium fraction of Bc/b-qnarks (see Ref. 0] for details). At 
T = 160 GeV, this fraction drops to as low as a few percent, but it is at least a factor of 
100 above the no-deconfinement scenario. To get a more realistic limit, we repeated the 
calculation, using only the initially-produced c quarks. This actually produces a much 
larger bound state fraction, since the c quark density only decreases with a volume 
factor rather than the exp(— rric/T) of chemical equilibrium. The last three rows in the 
Table show the corresponding B^. numbers at RHIC in this scenario. They depend quite 
strongly on the initial temperature, which determines the final charm density through the 
assumed isentropic expansion. 

We are in the process of refining these preliminary numbers. The kinetic equations will 
be followed numerically, using exact time evolution of the deconfined region and the effect 
of approach to chemical equilibrium for the charm quarks. It appears that the sensitivity 
to the parameters of the deconfined state will remain, making the observation of any 
BcS at RHIC both a "smoking gun" signal of deconfinement and a probe of the initial 
temperature of the system and the initial density of deconfined charm. 



Observable 


First Year 


Design Luminosity 


cc-pairs 


2.810^ 


6.5 10^ 


66-pairs 


1.210^ 


3.2 10^ 




1.610^ 


3.9 10^ 


T(ls) — > 


140 


3800 


Be J/iplu ^ fi^fi^lu 




(No Deconfined Phase) 


0.05-0.18 


1.5-4.9 


(QGP+cc in Chemical Equil.) 


18 


490 


(Only initial cc at To = 500 MeV) 


130 


3530 


(Only initial cc at To = 400 MeV) 


235 


6420 


(Only initial cc at To = 300 MeV) 


475 


12900 
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4.7. R. Vogt: Predictions for J/Psi Suppression 

The predictions are based on the model for J / ip absorption by nucleons, comovers, and 
quark-gluon plasma discussed in the recent review As in Ref. we take a nucleon 
absorption of 4.8 mb for color octet cc pairs and assume that the final-state charmonia can 
be broken up by mesonic comovers with an 0.67 mb cross section. The comover density 
is assumed to be proportional to the participant density in Au+Au collisions, tiauAu- 

The scale of the suppression is set by the transverse energy distribution of lepton pairs, 
shown in Fig. 32. In this calculation, the average transverse energy at a given impact 



parameter is the sum of hard and soft components and is approximately proportional 
to the number of collisions. Nuclear shadowing is included in the estimate of the hard 
component, see Ref. for details. The distribution is plotted as a function of -ErZ-E'T.max 
to be detector independent. Note that -Er.max ~ 1600 GeV for STAR and 700 GeV for 
PHENIX with an uncertainty of 10-20% depending on the shadowing model. 




Figure 32. Lepton pair Et distribution for RHIC. 



The resulting J ftp survival probabilities are shown in Fig. Figure ^(a) is a com- 
pletely hadronic scenario where the comover density is proportional to: nAuAm as for 
NA38 S-l-U; 2nAuAu5 as for NA50 Pb+Pb; and SwauAu, beyond reasonable expectations 



for hadronic matter. In Fig. |33|(b), plasma production is assumed in addition to comover 
absorption with density tt-auAu- Plasmas with 2 and 4 quark fiavors are studied. The 
energy density as a function of Et is determined from the hard and soft components to 
the average Et 0. When n/ = 4, only the Xc and ijj' are suppressed, but for = 2, 
direct J/ip suppression occurs when Et/Et. max ~ 0.2. Note that plasma is created in 
even the lowest Ets produced at RHIC so that Xc and ip' suppression begin immediately 
in both cases. 

The J/ip and Drell-Yan cross sections were calculated in perturbative QCD with and 
without nuclear shadowing. In the lepton pair mass range 4 < m < 9 GeV, Bajj^/ ajyy ~ 
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Figure 33. (a) J/ip suppression by hadrons only for comover densities proportional to 
'^AuAu (dashed), 2nAuAu (solid), and S^auAu (dot dashed), (b) J/ip suppression when 
quark-gluon plasma is produced for a plasma with 2 (dot dashed) and 4 (solid) quark 
flavors. The dashed curve is the same as in (a). 



190 — 280 depending on the shadowing parameterization. The Drell-Yan cross section is 
calculated at next-to- leading order. The J/tl) cross section has been calculated both in 
the color evaporation model and in the non-relativistic QCD approach. The two methods 
yield nearly identical results at RHIC when no nuclear shadowing is included. See Ref. 
for details of the cross section calculations and for more numbers. 
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5. EXOTIC POSSIBILITIES AT RHIC 

5.1. J. Schaffner: Ua{1) Restoration Signatures 

Strange Signals from the Chiral Phase Transition? [|l| 

Recently, there have been strong indications from lattice calculations, that the chiral 
f/^(l) symmetry is partially restored just above the critical temperature 0. It was found 
that the screening mass splitting between the S meson (the oq, a scalar-isovector meson) 
and the pion is reduced to about 5% or less. That means that the mass gap has changed 
by at least an order of magnitude. We are now going to discuss possible signals for 
RHIC associated with that effect. We explore the consequences with an SU(3) xSU(3) 
chiral Lagrangian. If the chiral SU(3) symmetry is effectively restored, then the masses 
of the pion and the sigma meson and the masses of the eta and the 6 meson are the 
same separately. If the chiral Ua{^) symmetry is effectively restored, then this mass gap 
vanishes and all four meson masses are the same: = wIq- ~ ms = for c ~ which 
gives two new observable effects for RHIC 0. 

Prediction 1: the number of r/s is enhanced by up to an factor four. The ao(980) meson 
has a width of F = 50 — 100 MeV. Its mass decreases with temperature, as seen on the 
lattice. Already below Tc the decay oq — + vr must be blocked by phase space (see 
figure). Also the matrix elements A^(ao — ^ r^ns + vr) = 4Acr in the chiral SU(3) phase 
and M.{aQ — >• 77^ + vr) = 2c ~ in the chiral Ua{X) phase. Hence, the inelastic channels 
are closed. The elastic channels are still large as they are proportional to the coupling 
constant A only. The ao can not decay to two pions which is forbidden by isospin. Also 
the decay to two kaons is heavily suppressed as the oq is actually lighter than one kaon 
alone. If the expansion from the chiral phase at Tc until the freeze-out temperature is fast 
enough (faster than the lifetime of the which is about 2-4 fm), then the numbers of 

~° is three times the numbers of the pions due to isospin counting. As the Oq decays 
mainly to r^s and pions, the numbers of the 77s can be increased by a factor up to four. 

Prediction 2: The scalar k {Kq) appears in the Kn mass spectrum. Recent studies of 
the ttK scattering data show the existence of a broad scalar, strange resonance, the fi;(900) 
[|]. The mass of this particle couples also to the Ua{X) anomaly, i.e. the mass splitting 
to the kaon is proportional to Amx oc c ■ a. The matrix element for the strong decay of 
the K to a pion and a kaon decreases in the chiral Ua{^) phase from values around F .8 
GeV to F < .2 GeV. Hence, the barely visible broad resonance gets a much smaller width 
in the chiral Ua{^) phase. As the mass of the k approaches that of the kaon, the strong 
decay k K + t{ is blocked by phase space already below Tc. Now, if the system freezes 
out around Tc, there appears a cusp structure in the Kir invariant mass spectrum between 
niK = ?^7r = -64 GeV and = .9 GeV due to the opening of the decay channel of the 
K. The cusp should be pronounced as the width of the k has substantially decreased. 

Prediction 3: strange clusters of two baryons are formed at RHIC. The status quo about 
the baryon-baryon interaction is, that there is some (very limited) knowledge about the 
NK (AS, AS, AA) interactions from A (S, H, AA) hypernuclei and AA, AE scattering 
and E~ atoms. There is absolutely no experimental information about the AS, AS, EE, 
ES, SS interaction at all. As hyperons are produced abundantly in relativistic heavy-ion 
collisions, this opens a new domain of strong interaction physics for RHIC. General SU(3) 
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flavor symmetry predicts that there are (quasi) bound states in e.g. S+p, S^E^, S^S^. 
The new Nijmegen soft-core potential fitted to scattering data indeed predicts bound 
states of EE, ES, and SS |^. These exotica are detectable by their weak two-body 
(nonmesonic) decay to final states like proton, A, and S~. Detectable candidates are e.g. 
lj:+p)b-^p + p, {E^p)b ^ A+p, (S°A);, E- +p or A + A, and (H-H")^ A + E'. The 
branching ratio is about (5 — 30)% depending on the channel and the binding energy [Q. 
The decay length is about cr = 1 — 5 cm, the more bound the candidate, the shorter the 
lifetime. The production rates have been estimated to be around 0.02 to 0.001 per single 
event using RQMD2.4 and wavefunction coalescence 
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5.2. D.H. Rischke: Parity Violation and Color Superconductors 

We give a pedagogical discussion of how, for massless quarks at very high densities, the 
formation of a spin- zero color superconducting condensate spontaneously breaks both the 
axial U{1) symmetry and parity |jl|. This observation is implicit in the seminal work of 
Baihn and Love, is noted by Afford, Rajagopal, and Wilczek, and is explicitly discussed 
by Evans, Hsu, and Schwetz 0. 

For simplicity, consider two degenerate flavors of quarks, and assume that a quark- 
quark condensate forms in the color-antitriplet channel |1|^- For massless quarks, two 
of the four possible condensates with total spin J = are |l|] 

{<Pi) = ^""'^ e/. ( C 75 ql ) and (0^) = e'^"^ ej, (qfciql), (18) 

where a,b,c = 1, 2, 3 are SU{3)c color indices, f,g = l,2 are SU (2) j flavor indices, and C 
is the charge conjugation matrix, (pl ^ are antitriplets under S'?7(3)c gauge transformations 
and singlets under SU (2) j rotations |jl|. The condensate 0" is even under parity, = 0"*", 
while 02 is odd, = 0^. There are two other condensates |1|, but they do not change 
our qualitative arguments about parity violation, and so we omit them. 

In the limit where mass and instanton-induced terms can be neglected, the effective 
Lagrangian for color superconductivity is 

>Co = \d^h\' + \dM' + ^ {\<Pif + I02|' - Ivl'^y , (19) 

where |0p = l]a(0")*0"- When mass and instanton effects are neglected, the Lagrangian 
is symmetric under axial U{1) transformations, which rotate 0" and 02 into each other. 
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Therefore, there is only one quartic couphng, A. The Lagrangian ([l9| ) generates nonzero 
vacuum expectation values for the 0"'s, which can be written as 

{(f)1) = cose , (0^) = v"" sine . (20) 

Condensation picks out a given direction in color space for f", and a given value for 
e. f " 7^ breaks the SU{3)c color symmetry, which produces color superconductivity. 
e ^ breaks the axial U{1) symmetry. Further, whenever e ^ 0, there is a nonzero 
= 0^ condensate (02) i ^^is represents the spontaneous breaking of parity (relative to 
the external vacuum). 

This breaking of parity is actually familiar from the spontaneous breaking of chiral sym- 
metry. Consider two flavors of massless quarks; the effective potential is 0(4)-symmetric, 
involving the = 0+ a- and = 0~ vr-meson fields. For massless quarks, it is as likely 
for a parity-odd pion condensate to form as it is for a parity-even cr-meson condensate. 
This does not happen in nature, because nonzero quark masses break chiral symmetry 
explicitly, and thus favor a 0"*" condensate. 

Similarly, it is important to add to the effective Lagrangian (|Ty|) terms which explicitly 
break the axial U{1) symmetry: 

C' = -c{\<p,\'-\<P2\')+m'\<P2f . (21) 

As shown by Berges and Rajagopal the first term is due to instantons, with c pro- 
portional to the instanton density. Instantons are attractive in the = O"*" channel, and 
repulsive in the = 0~ channel, so c is positive. 

In the second term, each power of the current quark mass rriq is accompanied by one 
power of 02- Since 02 itself is not gauge invariant, the simplest gauge- invariant term is 
mg|02p so m ~ niq. Thus, the pseudo-Goldstone boson for the axial U{1) symmetry is 
extremely light, m ~ 10 MeV, taking rriq to be the up or down quark mass and assuming 
the constant of proportionality between m and to be of order 1. This is in contrast 
to the explicit breaking of chiral symmetry, where the corresponding term is linear in 
the quark mass. The pseudo-Goldstone bosons are the pions which are relatively heavy, 
~ 140 MeV ~ ^/m^- 

Both instanton and mass terms act to favor the formation of the 0^ condensate 0i 
over that of the 0~ condensate 02- Consider, however, the limit of very high densities. 
When the quark chemical potential fi oo, the instanton density and so c vanish like 
~ (for two flavors). The real question is whether at some density the current quark 

mass is negligible compared to the scale of the condensate. If this happens, we reach an 
"instanton-free" region in which quarks are effectively massless, C can be neglected, and 
parity is spontaneously broken. 

Because mass terms are always present, the true thermodynamic ground state is always 
the parity-even 0"*" condensate, i.e., ^ = 0. There is, however, a finite probability for the 
system to condense in a parity-odd state, i.e., ^ 7^ 0. The size and lifetime of this state 
is set by the mass of the pseudo-Goldstone bosons. For chiral symmetry breaking, the 
characteristic scale is l/m,r ~ 1-4 fm. This is small compared to the time and length scales 
of a heavy-ion collision, so that parity-odd fluctuations average to zero. On the other hand, 
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the region in space-time over which a parity-odd color superconducting condensate forms 
is large, 1/m ~ 20 fm. If the collision time is shorter than this time scale, there is a finite 
probability that the system decays in a parity-odd state. We therefore propose to trigger 
on phase-space regions where nuclear matter is cool and dense, in order to observe the 
formation of parity-odd color-superconducting condensates on an event-by-event basis. A 
possible global parity-odd observable was discussed in 0. 
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5.3. D. Kharzeev: CP Violation in Au+Au? 

It has been proposed by R. Pisarski, M. Tytgat and myself [|I],0 that the discrete 
symmetries of strong interactions - parity V and CV - can be spontaneously violated in 
the vicinity of the deconfining phase transition. This would lead to a variety of dramatic 
effects which can be observed experimentally. What follows below is a brief and elementary 
introduction to the ideas of |[T|J^. I refer the reader to these papers for all details and 
further references. 

V and CV violation in strong interactions was never observed. However, in QCD, these 
symmetries cannot be taken for granted, and V and CV invariance should be regarded 
as a property of the ground state of the theory - the vacuum. When QCD vacuum is 
excited in a high energy nuclear collision, the properties of this excited vacuum state 
under V and CV transformations are in general not pre-determined, and will be defined 
by its structure. With these considerations in mind, let us recall briefly what is known at 
present about the V and CV symmetries of QCD. 

Classical equations of motion of QCD are known to possess topologically non-trivial 
solutions - four- dimensional configurations of the gluon field which are characterized by 
different values of the "winding number" n. In this situation, the true vacuum of the 
theory should be represented as a linear superposition of states with different n; this is 
analogous to the structure of the ground state in a crystal (Bloch wave function): 

\e)=Y.e''-\n), (22) 

n 

where 6 is called angle" , and n in terms of the gluon fields is defined as n = / df"x Q{x), 
with the topological charge density Q{x) = g^/{32TT^)tr{G^^G^''). Once an expectation 
value of a local observable is computed via the path integral, one has to include configu- 
rations of different winding numbers n with the weight given in (|2^). This procedure is 
equivalent to adding to the QCD Lagrangian a term 

^(^ = ^Q = ^^ HG,uGn- (23) 

Since in terms of the color electric, and color magnetic, fields, Gfj_yG^^ ~ ■ -B, it is 
easy to see that the "^-term" (^) explicitly violates V and CV invariances if 7^ 0. V and 
CV conservation in QCD thus is not guaranteed a priori] however current experimental 
constraints on 9 are very stringent, 9 < 10~^, and this constitutes so-called "strong CV 
problem" . 

In terms of an effective theory of Goldstone bosons, the 6'-term (^) manifests itself in 
the structure of the potential for the chiral U{Nf) matrix. By a global chiral rotation, 
a constant U field can be rotated into a diagonal matrix. With Uij = exp{i(f)i) 6'^\ the 
effective potential can be written down as 

m) = /.'(-cE^^co.(0.) + ^(E0.-^)') , (24) 



where the last term is proportional to the topological susceptibility, a ~ / d'^x {Q{x)Q{0)) 
- the correlation function of topological charge at zero momentum, and mi are current 
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quark masses. At zero temperature, the second term in (|2^) dominates, and the ground 
state is trivial, with = 0. 

Based upon an analysis in the limit of a large number of colors, we suggested [0] that 
near the phase transition, a becomes much smaller than its value at zero temperature. In 
a mean field type of analysis, with the temperature of the deconfining transition, and 
t = {Td — T)/Td the reduced temperature, we found that c(T) ~ and a{T) ~ t, so 

that the relevant ratio, a/c, scales as a(T)/c{T) ~ t^/^. Once a/c gets small near T^, the 
first term in becomes comparable in magnitude to the second, and (|2^) admits non- 
trivial metastable vacuum solutions with J2i{4'i) 7^ 0- It is clear from (^) that non-zero 
values of (01 + 02 + 03) act like having a system with non-zero 6— angle, 6^0. These 
solutions therefore correspond to metastable domains of V— and CP— odd vacua. Once 
the system is heated above the deconfinement phase transition and cools down, it can be 
trapped in one of these metastable V— and CP— odd states. 

An obvious question is how this breaking of the discrete symmetry by a metastable 
bubble could be measured in nuclear collisions. As the bubble is odd under V and CV, the 
pions produced by its decay must also be in a state which is odd under these symmetries. 
For the collisions of nuclei with equal atomic number, as the initial state is even under V, 
the observation of a "P-odd final state must be due to parity violation. In we proposed 
measuring, on an event by event basis, a global variable which is odd under V: 

J = ^ (p+ X p_) ■ z , (25) 

where the sum includes all tt+tt^ pairs in a given event, and i is a fixed vector of unit norm. 
Various choices for this vector, as well as a general classification of V— odd observables, 
can be found in 0. Since the effective potential (0) is symmetric, J on the average should 
vanish when summed over many events, and the distribution in J should be symmetric 
with respect to zero. V— odd effects therefore will manifest themselves in a non-zero 
width of this distribution. 

It is easy to understand why P-odd bubbles induce non-zero values of J in a given event. 
In terms of the underlying gluonic fields, the P-odd bubbles arise from fiuctuations in the 
topological charge density, G^^G^^ . Consider the propagation of a quark anti-quark pair 
through a region in which G^^G^'^ E ■ B 0. HE and B both lie along the z direction, 
then a quark is bent one way, the anti-quark the other, so that {pg x pg) ■ z ^ 0, where 
Pq and Pq are the three-momenta of the quark and anti-quark, respectively. An estimate 
of the effect can be done by considering the topologically non-trivial solutions directly in 
terms of collective pion fields we find that the P— odd asymmetries can be relatively 
large, at least ~ 10^^. P and CP violation in nuclear collisions is therefore possible, and 
should be searched for. 
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6. M. Gyulassy: Concluding Remarks 

We have seen considerable variation in section one on the predictions of global ob- 
servables (Bass, Bleicher, Gassing, Drescher, Eskola, Wang). Much of the factor of two 
uncertainty in the initial entropy and transverse energy are inherent our ignorance of the 
initial conditions and especially its soft component. However, as emphasized by Eskola 
and Wang a large part of that uncertainty is also due to interesting controllable physics. 
That component, nuclear shadowing and anti-shadowing of the structure functions, can 
in fact be inferred from systematic p + A studies at RHIC (Wang, Lin) or in e + A when 
electron beams can eventually be aimed through RHIC (Kovchegov). The p + A experi- 
mental studies are mandatory prerequisites not only to reduce the uncertainties in the 
initial conditions at RHIC, but also from the fundamental goal of simply understanding 
the nuclear wavefunction on the light cone. At first, light ion data, e.g. Si + Au, will be 
essential to get a feeling for the size of those effects, but eventually a full complement of 
pp and pA data will be needed as at the AGS and RHIC to untangle the physics. 

The A,B dependence of asymmetric systems as well as multiplicity dependence is also 
key to understanding the deviations from phase space saturation. Fireball models as 
shown in section 2 suggest that hadro-chemical equilibration is approached in strangeness 
to a surprising degree, and most microscopic models cannot reproduce the observation 
without invoking novel concepts such as baryon junctions (Vance) and higher Casimir 
flux tubes (Sorge). However, the linear variation of {K)/{tt) with multiplicity in NA49 
for example suggests that nonequilibrium effects play an important role. Deviations from 
flavor equilibration predicted by Stachel will also be important to look for at RHIC. 
Those deviations can provide us a handle to learn about the physics of dense matter. 
They may be due to phase space coalescence (Zimanyi), small relevant transport cross 
sections, or interesting medium mass modifications near the chiral symmetry boundary 
(Rapp, Schaffner). Just as the interesting physics are deviation from OSCAR certified 
[0] transport theories, the deviations from fireball fits must be carefully scrutinized as 
emphasized by Rafelski. 

I cannot over emphasize enough the importance of gaining control over the initial con- 
ditions. As Schlei and Dumitru pointed out, the data can be fit with almost any transport 
or hydro model is we are allowed the freedom to dial in arbitrary initial condition. The 
goal of the RHIC program is however not to tabulate arbitrary initial conditions but to 
understand in detail the energy and A dependence of the initial as well as the final condi- 
tions. Only with knowledge of the initial conditions can we interpret collective observables 
such as directed, transverse, and azimuthal asymmetric flow and such as possible time 
delay via HBT as evidence for new physics. 

One collective probe emphasized by Dumitru from the hydrodynamic point of view is 
dE±/dy that decrease in general if work is done by the system. However, given the present 
uncertainties in the absolute height of " Mount RHIC" , I suggest that the following simpler 
relative measure of longitudinal collectivity: 

R{t) = dET/dri/dNch/dr] 

Since Bjorken's work, it is known that an ideal p = e/3 equations of state leads to 
R{oo)/R{to) ~ 1/2 due to pdV work associated with longitudinal Bjorken expansion. 
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Covariant transport theory as discussed by Zhang and Molnar predicts less but still sig- 
nificant longitudinal work R{oo)/R{tq) ~ 3/4, , due to finite size dissipation. 

One of the most interesting observations at the SPS in my opinion was in fact NO evi- 
dence for any longitudinal work. Both E± and N^h are observed to scale almost perfectly 
with the number of wounded nucleons. Models such as RQMD and UrQMD provide a 
possible answer to the missing work puzzle. Heavy resonances and strings in such models 
mimic the mixed phase in QCD where the speed of sound is anomalously small. The 
mixed phase is notoriously lazy, as Hagedorn taught us! Excitation energy is wasted 
making heavy resonances instead of converting into collective motion. At RHIC, much 
higher initial energies are expected to be produced via mini-jets than at SPS, well into 
the plasma phase. There the speed of sound is expected to approach 1 / ^/3 again, the 
freshly produced plasma can push its neighboring cells more efficiently down the beam 
pipe with a resulting in decreasing R. 

It is important to distinguish this type of collectivity from that observed associated 
with radial and elliptic fiow and discussed in section 3. Those may also be sensitive to 
interesting variations in the equation of state. However, the longitudinal work occurs 
mostly during the early times when the longitudinal gradients ~ 1/r are the largest. 
Transverse collectivity may develop over longer times because the transverse gradients 
are ~ ^/Ra- The system must spend its first three fm deep in the plasma phase to get 
longitudinal work going. At the SPS lazy mixed phase produces no significant longitudinal 
work. Note that the hadron systematics in section 2 indicate that AA at those energies 
cross the transition region but the simplest signal of longitudinal collectivity was not seen. 
So my prediction is that we will start to see R decrease for the first time at RHIC as a 
function of centrality and A. 

One of the most exciting new frontiers at RHIC is high p± nuclear physics (see Eskola, 
Wang, Vitev, Snigirev, Srivastava). Unlike at SPS, where much theoretical ambiguity in 
pQCD results from uncertainties associated with intrinsic k± smearing as discussed by 
Wang, at RHIC collider energies the pQCD power law tail is predicted to stick out clearly. 
This provides a calculable, reliable base or calibration point from which deviation can be 
used to extract the physics of dense matter. In Fig. 34 from [|l| the contrast between 
the now known SPS and the predicted RHIC domain is evident. At the SPS HIJING 
accidently fits the data, but variations of uncontrolled soft hadronization assumptions 
leads to extreme discrepancies. At RHIC on the other hand, that soft component has 
much smaller infiuence on the high p± > 3 GeV domain. Therefore it becomes possible 
to look for new physics such as the non-linear energy loss of gluons predicted by BDMSP. 
The magnitude of predicted jet quenching phenomena is illustrated in the figure. 

Heavy fiavor physics at RHIC will also be interesting as emphasized in section 4. The 
open questions left by NA50 can be clarified at RHIC energies which is well above the 
charm thresholds and the high energy densities occur already in much smaller A systems. 
The anomalous J/tp deficit has been predicted to occur by Satz already in Cu + Cu, and 
brand X comover alternative scenarios will be much more easily ruled out as emphasized by 
Vogt. The topic of open charm is also rich with speculation and possibilities as discussed 
by Lin and Thews. 

Finally, I comment on section 5. The predictions for truly exotic and novel dynamics 
associated with chiral restoration, P violation at high baryon densities and especially 
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Figure 34. Opening of the high window at RHIC from [|I| . a) The WA80 S + S data 
(triangles) and the WA98 Pb + Pb ^ 7r° data (dots) are compared to HIJING1.35 with 
soft p± kicks (full lines) and without pt kicks (dot-dashed curves). The later scale with 
the wounded projectile number times a^A times the invariant distribution calculated for 
pp. The parton model curve from X.N.Wang is labeled by 'Wang'. The filled squares show 
pp 7r+ data scaled by the (Glauber) number of binary collisions times aAA for both SS 
and PbPb. b) Significant jet quenching predicted for RHIC energies 0, is masked by soft 
physics at SPS energies in the HIJING model. 



possible CP violating domains add great excitement to exploratory experiments at RHIC. 
While the predictions are based on very bold and even further extrapolations from known 
physics than those discussed in the previous sections, those directions associated with 
subtle many-body quantum phenomena should also be pursued vigorously because the 
potential payoff is so high. 
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